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This thesis is the result of my work as research fellow at IASF-MI, Milan section of 
the Istituto di Astrofisica Spaziale e Fisica Cosmica, part of INAF, Istituto Nazionale 
di Astrofisica. This work started in January 2006 in the context of the PhD school 
program in Astrophysics held at the Physics Department of Universita' degli Studi di 
Milano under the supervision of Aniello Mennella. 

The main topic of my work is the software modelling of the Low Frequency Instrument 
(LFI) radiometers. The LFI is one of the two instruments on-board the European Space 
Agency Planck Mission for high precision measurements of the anisotropies of the Cosmic 
Microwave Background (CMB). 

I was also selected to participate at the International Doctorate in Antiparticles 
Physics, IDAPP, IDAPP is funded by the Italian Ministry of University and Research 
(MIUR) and coordinated by Giovanni Fiorentini (Universita' di Ferrara) with the ob- 
jective of supporting the growing collaboration between the Astrophysics and Particles 
Physics communities. It is an international program in collaboration with the Paris PhD 
school, involving Paris VI, VII and XI Universities, leading to a double French-Italian 
doctoral degree title. 

My work was performed with the co-tutoring of Jean-Michel Lamarre, Instrument 
Scientist of the High Frequency Instrument (HFI) , the bolometric instrument on-board 
Planck. Thanks to this collaboration I had the opportunity to work with the HFI team 
for four months at the Paris Observatory, so that the focus of my activity was broadened 
and included the study of cross-correlation between HFI and LFI data. Planck is the first 
CMB mission to have on-board the same satellite very different detection technologies, 
which is a key element for controlling systematic effects and improve measurements 
quality. 

The thesis is organised in four chapters: 

1. a short introduction focused on state-of-art CMB phenomenology 

2. a chapter about the Planck mission mainly focused on the LFI instrument 

3. a chapter about software modelling of the LFI radiometers which includes a de- 
tailed description of the LFI instrument. Here I also discuss the experimental data 
available from the measurements campaigns on radiometer components, the model 
implementation and its validation against the frequency response measurements 

4. a chapter about the satellite thermal environment, with particular reference to 
the stage cooled at 4K, which is of key importance for both instruments. In 
this chapter I show the result of the analysis of the propagation of temperature 
fluctuations through the HFI. 

5. a chapter about cross-correlation of HFI and LFI data. In this chapter I describe 
the implemention of data analysis sessions in the KST data visualization software 
with the purpose of simplifying and standardsing the cross-correlation analysis 



1. Cosmology 



In this chapter I will give a short overview of the current knowledge on cosmology in order 
to understand the fundamental role of Cosmic Microwave Background measurements. I 
prefer not to follow the historical development, as usual, but focus just on the current 
state-of-art. 

1.1. Standard Model of Cosmology 

In the last twenty years, the collection of a huge amount of observational data has greatly 
contributed to test different theoretical models of the birth and evolution of the Universe 
and led to the definition of a generally accepted Standard Model of Cosmology. 

According to this model the Universe began about 14 billion years ago when it started 
expanding from a hot and dense state. It was highly homogeneous with fluctuations of 
the energy density which have subsequently grown by gravitational instability to form 
the cosmic structures (galaxies, galaxy clusters and superclusters) observed today. 

Observations picture a Universe uniform at large scale, this supports the Cosmological 
Principle which states that the Universe is spatially homogeneous and isotropic at large 
scale. Under this assumption it is possible to build a spherically symmetric metric model 
for the spacetime metric, called Friedmann-Lemaitre-Robertson- Walker (FLRW) metric, 

M- 

dr 2 

ds 2 = -dt 2 + a 2 (t)[- + r 2 {d(j) 2 + sin 2 6d(f))} (1.1) 

J. KT 

where: 

• r is the comoving distance, i.e. the distance in a reference frame expanding together 
with the Universe, 

• a(t) is the cosmic scale factor, which describes the expansion of the Universe 
and represents the rate at which two points of fixed coordinates (ri,9%,(f)i) and 
(r2, #2) <^2 ) increase their mutual distance with time. a(t) is adimensional and usu- 
ally normalised to be equal to 1 today (to), therefore a generic distance r evolves 
following this equation: 

r = a{t)*r(t ) (1.2) 

• k is the curvature parameter which can have only three possible values: k = for 
a spatially flat, Euclidean Universe, k = +1 for a closed Universe with positive 
spherical curvature, and k = —1 for an open Universe with negative hyperbolic 
curvature. 
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In order to study the dynamics of the homogeneous expanding Universe in the context 
of General Relativity [37] , it is necessary to solve the Einstein Field Equation |14j , whose 
modern formulation is: 



= SttGT^ (1.3) 

few details about this equation: 

• it assumes units where c, the speed of light, is equal to 1 

• G^, the Einstein tensor, expresses spacetime curvature 

• G is the universal gravitational constant 

• T^, the stress energy tensor, includes also the contribution of Dark Energy or 
vacuum energy (|38j): a hypothetical exotic form of energy that permeates the 
space and tends to increase the rate of expansion of the Universe; it is the favoured 
hypothesis for explaining the acceleration rate of the Universe expansion |31j and 

The simplest model of vacuum energy is the cosmological constant (|H1), 



whose stress energy tensor is written in equation 1.4 



Therefore Einstein field equation, 1.3, states that the curvature of spacetime is due 
to the matter/energy content of the Universe, which includes Matter, Radiation and 
vacuum energy. 



FLRW metric (Equation 1.1 ) has an analytic solution to into Einstein's field equations 
1.3 by assuming the energy momentum tensor T^ u isotropic and homogeneous. The 



results are called Friedmann equations: 



,d(i) 2 8?rG kc 2 
'a® ~ ~3~ P ~a 
a{t) AttG 



u^ a = — (L5) 



a(t) 



(p + 3p) (1.6) 



where p is the global energy density of the Universe, which includes contribution of 
Matter, Radiation and Vacuum energy, k is the curvature parameter defined above and 
G the universal gravitational constant. 

Cosmological constant, or vacuum energy, is the simplest model of Dark Energy, a 
hypothetical exotic form of energy that permeates the space and tends to increase the 
rate of expansion of the Universe; it is the favoured hypothesis for explaining the recent 
observations that the Universe appears to be expanding at an accelerating rate. 

Applied to a fluid with a given equation of state, the Friedmann equations give the 
time evolution and geometry of the universe as a function of the fluid density. 
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1.1.1. Hubble parameter 

Following FLRW metric, equation two points at a distance d = a(t)r will move 
apart with a velocity v = ar = Hd. This result was obtained by Hubble ( [19J ) thanks 
to measurements of the recession velocity of galaxies, and H is Hubble constant]^} 

H = 70.1 ± l.Skm/s/Mpc (1.7) 



Based on the first equation of Friedmann (1.6) it is possible to define the critical 
density 

3H 2 , . 

*=5R? (L8) 
such that density values p above, below or equal to p c refer to closed, open or flat 
Universe. 

p c is also used to normalise the densities in order to obtain relative parameters: 

• Matter: Vt m = — 

Pc 

• Radiation: f2 r = — 

Pc 

• Cosmological constant: £l\ = — 

Pc 

• Total: $7 = — = VL m + Cl r + £l\ 

Pc 

therefore = 1 means flat, > 1 closed and < 1 open Universe. 

Observations |27j strongly favour a flat Universe, with a global density near to the 
critical density. 

The equation of state, i.e. the relationship between energy density and pressure, can 
be written, for c = 1, simply p = wp. Once it is specified, the two Friedmann equations 
can be solved for the scale factor a(t) and give solutions for a Universe dominated by 
different type of energies: 

• Matter domination, w = -, p oc a(t) -4 , a{t) oc i 2 / 3 

• Radiation domination, w = 0, p oc a(t) -3 , a(t) oc i^ 1 / 2 ) 

• Cosmological constant domination, w = —1, p constant, a(t) oc e^ Ht ' 

Extrapolating the components time dependence back in time, it is clear that the early 
Universe was Radiation dominated; it followed a Matter domination epoch and, thanks 
to measurements of the present value of f^A [|^7j] , it is possible to state that just in the 
present epoch the cosmological constant and the Matter are at the same energy density 
level. 



1 Measurement from WMAP5, Baryon Acoustic Oscillations and SuperNovae, considering Lambda CDM 
model 
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1.1.2. Hot Big Bang 

Friedmann's equations and galaxies recession, propagated back in time, suggest an initial 
hot and dense state, named Big Bang [32J , whose physics is still beyond our understand- 
ing, which was rapidly expanding and cooling. 

Galaxies recession velocity measurements are a strong evidence of the Big Bang, in 
this section I will review the other main pieces of evidence of the Hot Big Bang theory: 

• Few minutes after the Big Bang, after the primordial plasma cooled and formed 
protons and neutrons, light elements started forming by nuclear fusion until the 
Universe became too cold to allow nuclear reactions. The primordial abundances 
of light elements (Helium-4, Helium-3, deuterium and lithium-7) can be computed 
from mathematical models as ratios to the amount of hydrogen, H. Mass ratios pre- 
dicted with this method agree with the abundances measured in galactic surveys: 
hydrogen ~ 76%, helium ~ 24%, see [9]. 

• The third evidence of a Hot Big Bang is galactic evolution and distribution: pop- 
ulation of stars have been evolving so that distant galaxies, observed as they were 
in the early Universe, appear very different from nearby galaxies, observed in a 
more recent state. Moreover, Big Bang simulations match well the observed star 
formation, galaxy and quasar distribution. 

• The evidence that, at least historically, had been the most important is the observa- 
tion of the Cosmic Microwave Background, a nearly isotropic microwave radiation 
at few Kelvins which was started propagating when the Universe was very hot. 
Its photons started propagating 379000 years after the Big Bang at a temperature 
of few thousands Kelvins. After ~ 14Gyr propagating in the expanding Universe 
these photons have lost most of their energy due to spacetime dilation and are now 
detectable as microwave black body at ~ 2.7i^ 

1.1.3. Inflation 

The standard Big Bang (BB) model however, cannot explain some features of the present 
Universe, in particular: 

• Flatness: Big Bang expansion cannot constrain the Universe to be so near to 
flatness at the level of few percent unless very special initial conditions are chosen 

• Homogeneity and isotropy: CMB photons, since the birth of the Universe, had 
time just to cover a small part of the Universe (~ 10 14 m), about 1° in the sky 
today, leading to just small part of causally connected regions. This contrasts 
with the fact that CMB is highly homogeneous at all scales. 

• Structure formation: the standard BB model provide no explanation of the source 
of the inhomogeneities which lead to structure formation 
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All these issues have been elegantly explained by inflation |17j . the theory which 
assumes the early Universe passed through a phase of exponential expansion driven by 
a negative-pressure vacuum scalar field, the inflaton. 



The flatness issue is solved by rewriting the Friedmann equation (1.6), as: 



<"" 1 - 1 '"» 2 = w ^ 

the right hand side of the equation is constant, therefore during inflation, the pa 2 
term increases extremely rapidly as the scale factor a grows exponentially, therefore the 
term |f2 -1 — 1| must decrease with time. This means that whatever is the initial value 
of the Universe energy density f2, after inflation it is forced to be very close to 1, an 
almost perfectly flat Universe. Subsequent evolution of the Universe will cause the value 
to grow, bringing it to the currently observed value of about 1.01. 

Homogeneity and isotropy are explained by the fact that the Universe was causally 
connected and therefore in thermal equilibrium just before inflation. 

Finally the inhomogeneities needed for structure formation are due to quantum fluc- 
tuations in the primordial plasma which are amplified to the seeds of the current astro- 
physical structures. 



1.1.4. Lambda Cold Dark Matter Model 

Lambda Cold Dark Matter model, or ACDM, is the Universe cosmological model which 
currently explains better the observations: it includes the Friedmann metric, the Hot 
Big Bang model and the inflation as explained in the previous sections. 
Its most important features are: 

• the Cold Dark Matter (CDM) model, where the dark matter is predicted to be cold 
(i.e. having a non-relativistic velocity at the epoch of radiation-matter equality), 
non-baryonic, dissipationless (cannot cool by radiating photons) and collisionless 
(i.e. interacting only through gravity). 

• the cosmological constant, A, a dark energy term that allows for the current ac- 
celerating expansion of the Universe 

The ACDM model is characterised by six basic parameters, which needs to be mea- 
sured through observation of the CMB, SuperNovae and Galaxies Sky Surveys. 



1.1.5. Content of the Universe 

Figure |1.1| shows the current energy content of the Universe derived from observations: 
atoms are only a small part, less than 5%, while Dark Matter and Dark Energy dominate 
the picture. 
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Figure 1.1.: Universe energy composition today 



1.2. Cosmic Microwave Background 

CMB is a microwave radiation which pervades all the Universe, it has a thermal black- 
body spectrum at the temperature of ~ 2.7 K, therefore its spectrum peaks at about 160 
GHz. 



1.2.1. Origin 

In the first 300, 000 years after the Big Bang, thermodynamic equilibrium between 
matter and radiation was maintained by Thomson scattering between free electrons and 
photons; after this period the temperature was low enough (~ 3000-fT) to allow for the 
combination of electrons and protons into neutral hydrogen: the reduced density of free 
electrons made the matter transparent to radiation, which started to propagate freely. 

The time when this event happened is called recombination or decoupling (about 
390,000 years after the Big Bang); the photons we receive today have travelled freely 
from a surface, called last scattering surface, to the Earth. They travelled into a an 
expanding spacetime and were redshifted from the original 3000-fT to ~ 2.7 K, and their 
temperature will continue to drop as the Universe expands. 



1.2.2. Anisotropies and power spectrum 



As described in section [LL3J inflation amplified quantum fluctuations and generated the 
inhomogeneities in the primordial plasma that subsequently seeded structure formation. 
The best measurements of CMB anisotropies until now was given by the satellite 



WMAP after 5 years of observations, figure 1.2 it shows a very detailed and complex 
picture of the inhomogeneities of the last scattering surface, their amplitude is: 



AT 



10 



-5 



;i.io) 



where T is CMB temperature and AT anisotropies temperature amplitude. 

The properties of CMB anisotropies are usually represented by expanding it into 
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Figure 1.2.: Temperature map of CMB anisotropics measured by WMAP after 5 years 
of mission 



spherical harmonics: 

AT 

— (M) = X)W/m(M) (1-11) 

where Y£ m (9, </>) is the spherical harmonic function of degree £ and order m, £ ~ n/8 is 
inversely proportional to the angular scale 9, and ag_ m are the multipole moments which 
relate the amplitude of the temperature fluctuations to the angular scales. Due to the 
cosmological principle, temperature anisotropics do not have a preferred direction in the 
sky therefore a^ m can be averaged over m to obtain: 

Ci = (\a em \ 2 ) (1.12) 

The CMB angular power spectrum is constituted by the coefficients Ci and it is the 
most important tool for comparing theoretical predictions with CMB measurements. 



Primary anisotropies Primary anisotropics originated during or before recombination 
and were generated by quantum fluctuations amplified by inflation. The inhomogeneities 
present in the Universe after inflation produce anisotropies through three different effects, 



see figure 1.3 



Acoustic oscillations: gravity and radiation pressure triggered sound waves that 
alternatively compressed and rarefied regions of the primordial plasma. After the 
Universe had cooled enough to allow the formation of neutral atoms, the pattern 
of density fluctuations caused by the sound waves was frozen into the CMB. The 
first peak in the CMB angular power spectrum (the so-called Doppler peak, at 
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Figure 1.3.: Cosmic Microwave Background power spectrum, main features 



about 1°) is therefore due to a wave that had a density maximum just at the time 
of last scattering; the secondary peaks at higher multipoles (lower angular scale, 
so-colled Doppler foothills) are high harmonics of the principal oscillations. The 
effect is directly proportional to the density fluctuations Ap: 

^c*^ (1.13) 
T p v ' 

• Gravitational perturbations: Photons coming from high-density regions un- 
dergo a gravitational redshift, this phenomenon is called Sachs- Wolfe effect (SW). 
Its effect impacts angular scales larger than the horizon at the last scattering 6 > 2° 
and I < 90° . The global effect of these anisotropies is 

^ = _^? (1 . 14) 

T 3 v ' 

where A<3? is the gravitational potential. 

• Silk damping: recombination is not an instantaneous process but takes a finite 
time, during this time the photons had time to diffuse and smooth out anisotropies 
at multipoles higher than ~ 1000, corresponding to angular scales 6 < 10' 

Secondary anisotropies Secondary anisotropies are due to interactions that take place 
between the last scattering surface and the Earth, the most important are: 

• gravitational lensing: CMB photons deviation by gravitational attraction of mas- 
sive objects like galaxies or clusters of galaxies 
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Sunyaev-Zeldovic (SZ) effect: Compton scattering of CMB photons by non-relativistic 
electron gas present in galaxy clusters 



Dipole anisotropy The Solar System is travelling at about 360Km/s with respect to 
its Local Group, this is evident and was first measured by the Doppler effect on the 
CMB which gives a strong signal, at about 3mK, superimposed to the primordial and 
secondary anisotropics. 

The dipole anisotropy signal was measured with high precision by COBE DMR (|25j) 
at the level of ~ 0.4%; thanks to its stability, strenght, availabilty through the full sky 
and the same spectrum as the CMB, the dipole is the natural calibrator for anisotropies 
experiments. 

Thanks to the calibration of the dipole, Planck is expected to have between 1 and 2% 
photometric gain uncertainties for observation lenghts shorter than a day, see [7] 



1.2.3. CMB Polarisation 

Thomson scattering of the CMB photons during recombination generated linearly po- 
larised photons due to the presence of quadrupole anisotropy. 

Only photons which last scattered in an optically thin region could have possessed a 
quadrupole anisotropy, and this depends on the duration of the last scattering. For a 
standard thermal history, the expected level of polarization is 10%, therefore the level 
of E-modes is about one order of magnitude less than temperature signal. 

Moreover, polarisation data are complementary to temperature measurements, and 
could be useful in breaking parameter degeneracies and thus constraining cosmological 
parameters more accurately. 



1.2.4. Foregrounds 

Foregrounds is a general term for referring any emission that confuse the primordial 
CMB signal after the time of last scattering. It contains extragalactic discrete point 
sources and galactic emissions due to synchrotron, free- free and dust. 

The best way to understand the impact of foregrounds in CMB measurements is fig- 
ure 1.4 where all 5 temperature maps from WMAP five years are shown. Galactic 
foreground emission is very strong and widespread, mostly at low frequencies. Sophisti- 
cated component separation techniques are used in order to disentangle primordial CMB 
from foregrounds. 

The main help for this task comes from the different spectral indexes of each foreground 
with respect to the CMB and to other foregrounds; the spectral index (3 is the exponent 
of the power law dependence for brightness temperature Tf, with respect to frequency v: 

T b = v-P (1.15) 

multifrequecy measurements are therefore mandatory for separating the different con- 
tributions. 
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Figure 1.4.: WMAP 5 years temperature maps at all frequencies after dipole removal 



Figure [1. 2. 4 displays the application of this concept to the Planck space mission, that 



will be presented in section [2| Planck has a huge frequency spectrum coverage and it 
will measure parts of the spectrum where the relative contribution of the components is 
very different, the extreme case is 857 GHz which is going to measure the distribution 
of dust emission. 
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Figure 1.5.: Synchrotron, dust, and free- free foregrounds spectra superimposed to 
Planck channels bandwidths and to the expected CMB signal at 1 degree 
scale; in log-log scale, the exponential dependence of foregrounds emission 
versus frequency is represented by straight lines 
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1.2.5. Precision Cosmology with CMB 



A precise measurement of the CMB power spectrum allows to validate the prediction 
of different cosmological models, as the KCDM model reviewed in section 1.1 A\ and to 
compute with high precision the related cosmological parameters, six for the AC DM 
model. 

The most simple cosmological parameter measurable with CMB is the total energy 
density parameter f2, it is strictly related to the position (I) of the first peak; analytically 
can be computed with the following equation: 



peak 



220fT 1/2 



(1.16) 



The more general approach, which is used to estimate most of the cosmological param- 
eters, is the comparison of the measured power spectrum with predictions given by the 
cosmological models. For example simulations of CMB power spectra for different value 



of f2 are shown in figure 1.6 a measurements with low uncertainty is able of constrain 
more efficiently the value of f2. 

Equivalent procedures can be used to study other cosmological parameters, but more 
often all the relevant parameters are evaluated simultaneously using Monte Carlo tech- 
niques. 




Figure 1.6.: CMB power spectrum simulation for different values of Q 
The accuracy of the power spectrum estimation can be analytically modelled consid- 
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ering instrument noise and angular resolution, see |26j : 



SCi _ 1/2 / 2 A9p ix ap ix 



C e ~ fs ^ V 21 + 1 [1 + N ptx C e W E ] (L17) 



where: 



• the term is the contribution due to Cosmic Variance, which is the ultimate 

21 + 1 

accuracy limit for CMB power spectrum and is related to the fact that the CMB 
field is just a single realization of a stochastic process and it is not expected to 
follow the average of the possible realization, this effect is strong on the largest 
scales. 

• J sky is the fraction of sky covered, which has to consider also eventual cuts of the 
galaxy plane 

• 6pi x is the angular resolution of the optics 

• cjpix is the noise per pixel 

• Wg is the beam window function 

• Nm X is the number of pixels 

• A is the surveyed area 

Figure 



1.7 



SCg 

shows the values of — — for I = 1500, the scientific requirement of the 

Ce 

Planck mission is a precision of 5%, which is represented by the grey area in the figure; 
the rectangle encloses the expected performance of the 100 GHz Planck HFI channel. 
These requirements are strong drivers in instrument design: ~ 10 arc-minute angular 
resolution at 100 GHz requires an antenna of ~ X/9 ~ 1.5m, [iK sensitivity requires an 
array of feed horns cryogenically cooled with at least 1 year of integration time. 

This estimation didn't take into account calibration and systematic errors; their im- 
pact should be < 1 — 2%a p i x in order not to dominate the power spectrum; calibration 
accuracy of the order of 1% is achivable using the CMB dipole signal, systematics due 
to foregrounds need good frequency coverage, as explained above. 
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Figure 1.7.: — — for I = 1500 as a function of the angular resolution 6 and of the average 

w 

noise per pixel a p i x 
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2. Planck Mission 



The European Space Agency Planck satellite is designed to fully extract the cosmolog- 
ical information contained in CMB anisotropy by setting angular resolution, spectral 
coverage and sensitivity such that the power spectrum reconstruction will be limited 
by unavoidable cosmic variance and astrophysical foregrounds. Planck will also pro- 
vide a precise measurement of the TE correlation and of the E-mode polarisation power 
spectrum and possibly offer a first B-mode detection. 

Planck objective is to provide full sky maps with angular resolution between 5' and 
33', spectral coverage between 30 and 857 GHz and sensitivity of ~ 10~ 6 . Planck 
performance and design derive, as discussed in section 1.2.5, from the objective of 5% 
precision on CMB temperature power spectrum at small scales (1500 < I < 2000) 



Primary mirror 
Baffle 

Secondary mirror 




Focal plane 



Figure 2.1.: Overview of the Planck Satellite 



The satellite, see figure 2.1, will orbit around the Lagrangian point L2 spinning at 1 



rpm around the Anti Sun axis; the telescope is pointed at almost 90° with respect to 
the spin axis, therefore Planck will map the sky in circles; each circle will be scanned for 
an hour, 60 times, then the satellite will be repointed by 2.5'. 

The scientific payload is placed in the focal plane of an off-axis aplanatic telescope 
with a primary of 1.9x1. 5m, and it is constituted by two instruments for CMB detection: 
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• HFI, High Frequency Instrument, is an array of 52 bolometers cooled to 0.1 K 
between 100 and 857 GHz, |23] 

• LFI, Low Frequency Instrument, is an array of 22 pseudo-correlation differential 
receivers based on HEMT technology cooled to 20 K and operating between 30 
and 70 GHz, @j 

The Planck cryogenic chain, needed for cooling the scientific instruments, is composed 
by a passive cooling system and three independent stages of active cooling driven by 



different coolers, details about Planck thermal design are in section 4.1 

2.1. HFI 

The High Frequency Instrument (HFI) is an array of 52 receivers in 6 frequency bands (in 
the 100-857 GHz range) based on bolometric technology and cooled to 0.1 K. Bolometers 
are composed by a conductive grid and a piece of semiconductor in thermal contact with 
a cold bath; when a photon hits the bolometer, it changes its temperature and therefore 
its resistance. Its resistance is monitored by the readout electronics which records a 
signal. 

Two kinds of bolometers are used in HFI: 

• Twenty spider-web bolometers in the 143-857 GHz range, absorbing radiation via 
a spider-web-like structure 

• Thirty-two polarisation-sensitive bolometers in the 100-353 GHz range, absorbing 
radiation via a pair of linear perpendicular grids. (Each grid absorbs only one linear 
polarisation.) This kind of bolometer, called PSB, allows to measure polarised 
radiation. 

The temperature of 0.1 K is achieved by three coolers: a hydrogen Sorption cooler, 
[5], (providing a stage at 18 K), a Joule-Thomson mechanical refrigerator (precooled by 
the 18 K cooler and providing 4 K) and an open- loop 3 He/4 He dilution refrigerator, 
which provides 0.1 K. 

The 4K stage is very important because it is a key element in the HFI cryogenic chain 
and also provides a stable reference load for the LFI pseudo-correlation radiometers, see 
section I 



Figure 2.2 shows the position of HFI horns on the focal plane; HFI naming conven- 
tion is based on the frequency and an ordinal number for spider web bolometers (see 
section 2.1), polarization sensitive bolometers instead are labelled with the same ordinal 



number and a or 6 for distinguishing the ortogonally polarization sensitive channels. 
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Figure 2.2.: Schematic of the HFI focal plane with horns labels; the two channels of PSB 
bolometers, square in the schematic, are indicated as a and b, e.g. 100-lb 
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2.2. LFI 



LFI is an array of 22 pseudo-correlation radiometers (discussed in the following section 



2.3.3) connected to 11 horns and capable of measuring CMB temperature and polariza- 



tion in three frequency bands centered at 30, 44 and 70 GHz. 




Figure 2.3.: Picture of the LFI instrument during integrated tests in Alcatel Alenia Space 
Milano 



LFI, see figure 2.3 is separated in two sections: 



the Front End Unit, cooled by the Sorption Cooler at ~ 20K, includes the horns 



and implements the pseudo-correlation strategy as detailed in section 2.3.3 



• the Back End Unit, located into the satellite service module and kept at about 
300 K, includes a further amplification stage, the detection system, conversion to 
digital and compression. 

The microwave signal is brought from the FEU to the BEU with waveguides about 
1.5m long; the three V-grooves are passive radiative cooling system which favour thermal 
decoupling between the FEU and the BEU by setting the temperature of the interface 
with the waveguides at 60, 80 and 140 K. 
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After detection by the Back End Module, sky and reference load signals are integrated 
and digitized by the Data Acquisition Electronics (DAE). The digital signal is then sent 
to the Radiometer Electronics Box Assembly (REBA) which downsamples, compresses 
the data and builds packets that are sent to the Earth thanks to the transmission an- 
tenna. 



Figure 2.4 shows the position on the focal plane of LFI horns. 





Figure 2.4.: Planck focal plane with LFI horns labels 



My work is mainly focused on LFI, therefore general description of the radiometric 
technology and radiometers frequency response will be presented in the next sections, 
while the details of each LFI component are presented in the chapter about the radiome- 
ter model, see chapter [3j 



2.3. Radiometric technology 

In this section I will give a brief introduction on the available types of radiometers in 
order to understand how the pseudo-correlation radiometer was conceived, then I will 
show the radiometer response analytical representation focusing on bandpass response. 
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Figure 2.5.: Schematic of a total power radiometer 



2.3.1. Total power radiometer 

Radiometers are coherent receivers which are able to measure the radiative power of 
electromagnetic waves. Radiometers operating in the microwave range, in their simplest 



form, called total power radiometer, are composed by, see figure 2.5 



• an antenna, usually a corrugated feed horn, to transfer the electromagnetic waves 
from vacuum or air to the transmission line reducing losses due to impedance 
mismatch 

• one or more amplifying stages, with a total gain in terms of power between 10 6 
and 10 9 , which corresponds to 60-90 dB 

• a passband filter, needed for select the interesting frequency band for the measure- 
ment 

• a detection stage (usually a square law detector, typically a diode) which trans- 
forms the Radio Frequency signal to Direct Current 

• an integrator, which integrates the DC signal in order to reduce fluctuations 
The sensitivity of a total power radiometer is given by: 

ot ^ y - T = — (2.1) 

where r is the integration time, i.e. the observation time, and Af is the bandwidth. 

2.3.2. Dicke switching radiometer 

Amplifiers gain fluctuations, usually called 1// because they increase at low frequency, 
are predominant on signal in this extreme applications; a possible strategy is using a 
stable thermal source as a reference and quickly switching between the antenna and the 
reference load. The Dicke radiometer is basically a total power radiometer with two 
additional features: an RF switch and a synchronous demodulator inserted after the 
square-law detector. 



The input of the radiometer (see Figure 2.6) is rapidly switched between the antenna 



temperature and the reference temperature. The switch frequency Fs is typically of 
the order of few GHz. The output of the square-law detector is multiplied by +1 or 1, 
depending on the position of the Dicke switch, before integration. 
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Figure 2.6.: Dicke switching radiometer schematic 



It is important that the switching frequency is quicker than the typical amplifiers gain 
fluctuations, in this case: 

Vout — Vsky ^ref — G(T s ky ^noise) G(T re j + T no i se ) — G(T s j~y T re j) (2-2) 

The paper [33] shows an analytical tratment of pseudo-correlation radiometers, show- 
ing that Dicke switching radiometers reduce strongly the impact of systematic effects. 
Moreover, if the T s k y and T re f are exactly the same, the radiometer is completely not 
sensitive to gain fluctuations. 



Different reference load temperature If the reference load temperature is different 
from the sky load temperature is still possible to fight 1/f noise with this strategy by 
multiplying the T re f temperature by a number, r, called gain modulation factor: 



r = (2.3) 

T re f 

where the overline indicates an average computed usually on a time period of one 
hour. 

In this case the diffenced output is: 

Vdiff = Vsky — fVref = G{T s k y + T no i se ) — rG(T re f + T no i se ) = G(T s ky — rT re f) (2.4) 

The main issue of the Dicke switching radiometer is that for half of the observing time 
the radiometer is receiving the signal of the reference load, and this highly impacts on 



its sensitivity, see equation 2.5 



a T ~ 2 Tsky t^l mse (2.5) 
VtAz/ 

Therefore a Dicke radiometer has a sensitivity two times the sensitivity of a Total 
power radiometer. 



2.3.3. Pseudo-correlation radiometer 

The evolution of the Dicke radiometer is the pseudo-correlation radiometers, in this type 
of radiometer, sky and reference are simultaneously observered, see for example Planck 
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LFI implementation: the upper half of figure [2771 the cryogenic Front End Model, imple- 
ments the pseudo-correlation, while a second amplification, filtering and detection stage 
is present in the wark Back End Module. The hybrid outputs the sum and difference 
which are amplified by two independent amplifiers; the second hybrid splits the amplified 
signals back into sky and reference. 

Both sky and reference streams were amplified by both amplifiers, and are strongly 
correlated; the correlation is due to the 1/f noise of the amplifiers, therefore by taking 
the difference, similarly to the case of the Dicke radiometer, it can be strongly reduced. 

The advantage of such a configuration is that there is no signal loss and the sensitivity 
is comparable to the sensitivity of the total power radiometer, in particular: 



~ V2 



(2.6) 



The problem of this configuration is that only the amplifying stage involved into the 
pseudo-correlation is stabilised; the 1/f noise of the next stages is not affected. For this 
reason between the amplifier and the second hybrid a phase switch phase lags at 4 KHz 
one of the arms of a radiometer. The effect of the switching is that the sky and reference 
outputs are exchanged continuously; based on the same principle of the Dicke switching 
radiometer, this technique is used to reduce the 1/f fluctuations of the Back End unit. 
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Figure 2.7.: Planck LFI radiometers, in the upper section the cold Front End Module 
which implements the pseudo-correlation, in the lower section the warm 
backend 
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2.3.4. Basic radiometric equations 

In this section I will review the basic analytic modelling of the response of a radiometer 
and then I will focus on its bandpass response. 

If we assume a top-hat frequency response in a band Av around the centre frequency, 
the voltage output from each of the 44 output diodes can be represented mathematically 
by this equation: 

Vout — G • Pin ~\~ V no i se (2-7) 

where Pi n represents the input RF power, G is a proportionality constant (the so-called 
photometric calibration constant) and Vnoise is a noise contribution from the RF ampli- 
fiers and back-end electronics. 

It is common to refer this noise contribution to the input, so that: 

Vout = G • (Pin + Pnoise) (2-8) 

where V no ise = G ' Pnoise ■ 

Pnoise can be converted into a temperature, the Noise Temperature T no i se : 

rp 1 noise q\ 

-'-noise — ~, I \^"") 
KbAv 

where ks is Boltzmann's constant, Av is the bandwidth, T no i se is the input load tem- 
perature that would cause the V no i se output signal in an ideal system, i.e. a system with 
no noise with V ou t = G • Pi n ). 

Noise temperature is directly related to radiometer sensitivity by the radiometer equa- 
tion: 

o-n = k T ^ + T ™ (2.10) 
VAi/ • r 

where: 

• o"tv is the root mean square of the white noise 

• K is a factor which is dependent on receiver architecture. For the LFI receivers 
we have K = -J2 (H3J). 

• Av is the bandwidth 

• t is the integration time 

Considering in-band behaviour 



If we now consider explicitly the frequency dependence, Eq. 2.7 reads: 



Vastly) = G(u) ■ P m (v) + V noise (v) (2.11) 
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where the function G(u), i.e. the gain versus frequency, is also referred to as the receiver 
bandpass. 

LFI square law detectors integrate the signal over the band, therefore each channel 
output signal is: 



i Aw 



Vout= / V out {u)du (2.12) 



Aw 
vq 2T 



The next equations show all the relations between integrated and frequency dependent 
quantities: 

i Aw 

f u o+ — 

Vnoise = / V n <rise(y)dv (2-13) 

/,, An 

JVQ- — 

Pin = / Piniy)dv (2.14) 

C + ZG{v).P m {v)dv 
G = -!LJJ (2.15) 

in 

G, the photometric calibration constant, is the conversion factor between a variation 
in radiometers voltage outputs we measure and the corresponding variation in the in- 
put temperature; LFI radiometers, as WMAP ones, provide, by design, only relative 
measurements of the sky signal. 

During Planck Spacecraft operations this calibration factor will be computed using 
the well known modulation of the Cosmic Microwave Background dipole signal among 
the whole sky and then used to produce temperature maps of the sky, see |1.2.2| 



As evident from G equation (2.15), calibration is tightly bound to the spectral in- 
dex, or frequency dependence, of the input source. Free-free, dust and synchrotron 
foregrounds, instead, due to their own spectral indexes, have a different coupling with 
radiometer bandpass and their photometric constant would be different. Therefore, us- 
ing the photometric constant computed using CMB signal introduces a second order 
systematic effect on these components. 

Knowledge about bandpasses can be used to reduce the effect of this calibration 
method by using an iterative process to separate different components. 

The receiver bandpass, G{y) can be determined by performing, for each frequency, 
two measurements at different input power levels and computing: 

<M - (2,6) 
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3. LFI radiometer model 



3.1. Introduction 

3.1.1. History of the model 

The first radiometer model was developed in 2002 by Paola Battaglia, [3], and subse- 
quently improved by Cristian Franceschet, |16j . 

The model has been developed using the software simulation platform Advanced De- 
sign System by Agilent Technologies, see ADS website, The Advanced RF model original 
aim was to simulate the impact of subcomponents non- idealities, e.g. insertion and re- 
turn losses, and of systematics, e.g. OrthoMode Transducer asymmetries and thermal 
fluctuations, on instrument performances. 

Advanced RF model started as a complete analytical simulator using components 
embedded into the Agilent Advanced Design system, configured with expected perfor- 
mance defined in Planck LFI specification documents [3J. Measurements of the true 
hardware performance were made then available by the manufacturers, therefore some 
of the analytical models were replaced with data coming from measurements. 

During the Planck Low Frequency Instrument Flight Test campaign held in Thales 
Alenia Space Italia (Vimodrone-Milan) all the measurements available channel by chan- 
nel were included in the model. This improved model was then used to compare sys- 
tematically simulations to measurements with two objectives: 

1. making a cross validation of the two methods available for estimating bandpasses: 
RCA frequency response measurements and simulations using bandpasses of single 
components 

2. obtaining, through simulations, data that were not available from measurements 

The model had several issues related to its implementation into ADS, mainly due to 
the fact that ADS is a closed environment, it is not possible to run batch simulations 
on all RCAs and data export must be performed manually. For these reasons it was a 
very long process to produce a set of simulations of all channels and it was possible to 
introduce errors. Moreover, being a closed source software, it was not possible to read 
the source code in order to understand exactly what was the implementation strategy 
of the embedded models further than the details provided in the documentation. 

While continuing the development on the Advanced LFI RF model, I started imple- 
menting the model from scratch on the electrical circuit simulator software QUCS [6j, 
website http : / / qucs . sf . net . 

The main objective of this work is to take full control of the model, in particular on 
two aspects: 
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1. automation - in ADS it is not possible to run many simulations together and export 
automatically the results, therefore the run operation is quite long and tedious 



2. sourcecode - Cristian Franceschet found a bug in ADS that was not solved, because 
Agilent wasn't interested in continuing the development of that piece of software. 
Having source code available allows the users to fix the software themselves if 
needed 

I named the new QUCS based model QIMP, which is the acronym for QUCS Integrator 
of Measured Performances, it became stable and superseded the Advanced LFI RF model 
in September 2008. 



3.1.2. Chapter structure 

The first section is a brief introduction to the QUCS software. 

The second section is devoted to review the model implementation component by com- 
ponent referring to the hardware setup and to the data available from single component 
test campaigns. 

The third part gathers all the results obtained by simulations. 



3.2. QUCS 

Quite Universal Circuit Simulator (QUCS), [6] is an open source electronics circuit sim- 
ulator software released under GPL license. It gives the possibility to set up a circuit 
with a graphical user interface using QT libraries and simulate the signal and noise be- 
haviour of the circuit. It is an alternative to well known Berkeley SPICE and Agilent 
ADS. Synopsis of QUCS internals is given in Figure [37T| 



QUCS is implemented in C++ and use extensively class inference in order to facilitate 
the implementation of new components. The main QUCS modules are: 

qucsator the command line circuit simulator, which reads a circuit description in a prede- 
fined ASCII format (named netlist) and outputs an ASCII format results file 

GUI Graphical User Interface is completely independent and makes it possible to draw 
a circuit using a library of devices or file defined components. The GUI can also 
automatically build the netlist, run qucsator and parse the results file, allowing 
the user to easily produce tables and plots. 
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Figure 3.1.: QUCS synopsis: the high level structure, at the bottom of the figure, 
shows the typical simulation process, in the top more insight into QUCS 
implementation 
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3.3. LFI Radiometer model implementation 



3.3.1. The Low Frequency Instrument - Overview 



For a complete introduction on radiometers please refer to section 2.3 

LFI is a radiometer array mounted in the focal plane of Planck satellite telescope, 
receiving Cosmic Microwave Background photons at 30, 44 and 70 GHz. The complete 
LFI is an array of 11 Radiometer Chain assemblies (RCA), 6 with a central frequency 



of 70 GHz, 3 of 44 GHz and 2 of 30 GHz. Each RCA is composed by (Fig. 3.2): 



• a Feed Horn and an OrthoMode Transducer (polarisation splitter) looking at the 
Sky signal 

• 2 short waveguides (just for 30 and 44 GHz, 70 GHz reference horn are directly 
connected to the FEM) connected to two reference loads at about 4 K 

• a cold (~ 20K) pseudo correlation stage (Front End Module) providing ~ 35dB 
amplification 

• 4 waveguides connecting the cold module to the warm backend 

• a second RF amplification, a diode and post-detection DC electronics stage at 300 
K (Back End Module) 




Figure 3.2.: Picture of a 30 GHz Radiometer Chain Assembly (RCA) mounted for testing 
in Milano 



LFI channels naming scheme is based on: 

• RCA number (18 to 28) 

• OMT arm (Main or Side) 
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Figure 3.3.: model schematic: FEM is modelled by 4 independent bandpasses, therefore 
each of the 44 LFI channels is modelled independently as the above model. 
Each component is represented by its measured behaviour as a function of 
frequency. Waveguides instead are simulated analytically given their dimen- 
sions, length and temperature gradient along their length. 

• radiometer number (which is for the radiometer connected to OMT Main arm, 
1 for Side) 

• detector number (0 or 1) 

for example LFI22M-01 refers to RCA 22, OMT Main arm, channel 01 among the 4 
RCA 22 channels (00,01,10,11). 

3.3.2. LFI Radiometer model - Overview 

QIMP (QUCS Integrator of Measured Performances) can simulate the in band behaviour 
(frequency response) of each RCA using simplified schematics and data of each hardware 
subcomponent made available by the manufacturers thanks to performance tests. 

Although in this thesis the bandpasses provided by the model are often referred to as 
"simulated" bandpasses, we want to stress the fact that they are indeed the product of 
a combination of real measurements performed on the various subsystems to estimate 
the complete RCA frequency response. 

The most important model result is the bandpass which is useful for building CMB 
polarisation maps and for foregrounds removal; but the model also simulates volt out- 
put as a function of sky input temperature, noise temperature, input power at each 
subcomponent, photometric gain and so on. 

The model is built to study the broadband behaviour of the instrument in the fre- 
quency domain, therefore it is not able to produce data streams as a function of a input 
temperature changing with time, but can simulate statically the output with different 
input temperatures. 

The complete model implementation is shown in figure pT3| the frequency dependence 
of all parameters is implicitly considered if not otherwise specified. 
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3.3.3. Sky Load 

Sky load is simply modelled as a noisy resistor with tunable temperature; it outputs an 
uncorrelated signal independent of frequency with a power spectral density of: 

Pskyioad = k B T [Watts/Hz] (3.1) 

where ks is Boltzmann's constant and T is the chosen input temperature; the skyload 
assumes perfect absorption, i.e. no emission from the horn is reflected back. 

3.3.4. Feed Horn and OrthoMode transducer 
Description 

The first part of each RCA is a corrugated feed horn which collects radiation coming from 
the telescope secondary mirror reducing reflection and absorption thanks to an optimised 
design. The OrthoMode Transducer receives the sky signal from the horn and splits the 
orthogonal polarisations, that propagate through two independent radiometers. 

Each LFI passive component frequency response has been characterised by tests per- 
formed at IFP-CNR (Istituto di Fisica del Plasma - Milano). 

Model 

The Feed Horn and OMT assembly in QIMP has been modelled as a single 2-ports 
S-parameters component: 

• The main contribution to reflection (Sll) is the OMT, due to the waveguides 
configuration needed for splitting the incoming wave polarisations, therefore the 
Feed Horn reflection losses can be neglected and the Sll tests on standalone OMT 
taken as representative of the behaviour of the FH-OMT assembly. 

• S22 of the assembly, i.e. the reflection losses on FEM side, was measured directly 
on Feed Horn and OMT together. 

• S21, which is the most important parameter, because influences how much sky 
signal can reach the radiometer, and includes both reflection and resistive losses; 
for some devices we have only S21 measurements and not S12; however, looking 
at devices where both measurements are available, they are very similar, therefore 
if S12 is not available, it is considered equal to S21. 

The Sll and S12 parameters characterise respectively the signal reflected back into 
the Horn by the OMT itself and the signal reflected by all the instrument components 
after the OMT which travels back towards the Feed Horn. This parameters are included 
in the model and can be used to perform a study on LFI emission (both in testing 
configuration, e.g. standing wave between the sky load and the OMT, or in flight, e.g. 
LFI interference on HFI), but for now they are not used. 
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Measurements extrapolation 

The nominal bandwidth of LFI receivers is 20% of the central frequency, e.g. for 30 
GHz the nominal bandwidth is 6 GHz. Moreover, the standard band of the WR28 
waveguides, used for testing is 26.5-40 GHz; for this reason 30 GHz OMTs were only 
measured between 26.5 and 40 GHz. 

However, the fact that the measured bandpass response of the RCAs is still high at 
26.5 GHz showed that the radiometer band extends below this limit. For this reason, 
one of QIMP objectives was to provide the best estimation for the band between: 

• 21.3 GHz, which is the low frequency cut of the 30 GHz waveguides 

• 40 GHz, where the RCA response is already very low due mainly to BEM lowpass 
filter 

OMT design was exactly the same for 30 and 44 GHz devices, it was scaled in order 
to work efficiently with different wavelengths, therefore it is possible to exploit this 
similarity and extrapolate the measurements on 30 GHz devices using as a reference the 
measurements on 44 GHz which were performed on a larger band, from 33 to 50 GHz. 

The extrapolation was done after normalising the frequency axis on the central fre- 



quency, see 3.4 so that all the curves has a central frequency of 1. 



3.3.5. Front End Module 
Description 

The Front End Module is the cryogenic amplification stage where the pseudocorrelation 



is implemented, see section 2.3 



It is constituted by (figure 3.5): 



• an hybrid which receives the incoming sky and reference signal and outputs the 
sum and the difference 

• a cascade of 4 Low Noise Amplifiers based on High Electron Mobility Transistors 
which amplifies the signal between 30 and 40 dB, i.e. up to a factor of le4 in power 
and 100 in voltage amplitude 

• phase switches: a system through which the wave can travel along two alternative 
pathways phase-lagged by 180°. The pathway is selected by two diodes that are 
alternatively biased at a frequency of 4 KHz; this fast switching is implemented to 
reduce y gain fluctuations of the Back End Module 

• a second hybrid which splits the signal back into sky and reference components 

Originally the Advanced ADS LFI RF model implemented the pseudocorrelation, 
therefore hybrid, phase switches and LNAs where modelled as distinct components. 
Ideally this configuration is the most realistic and allows to compute also isolation, i.e. 
the level of separation of sky and reference channels of the same radiometer. 
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Figure 3.4.: Extrapolation of the 30 GHz measurements performed by comparison with 
44 GHz measurements normalised on their own central frequency. The ex- 
trapolation was necessary between 21.3 GHz and 26.5 GHz 



However, FEM subcomponents are tightly integrated, and we cannot make the as- 
sumption that the behaviour of a sub component tested alone is the same when in- 
tegrated into the FEM; the effect of connection, matching and coupling is strong and 
FEM designers and builders judged that the tests on the complete FEM performed after 
integration are more reliable than using sub components measurements. 

A test on the complete FEM means measuring gain and noise temperature as a func- 
tion of frequency of each of the four output ports when a signal is injected into the sky 
input. 

Therefore, in QIMP, the FEM structure is simplified and consists in 4 independent 
channels, neglecting pseudocorrelation. This simplification doesn't impact on QIMP ca- 
pability to model bandpass response because frequency response is a property intrinsic 
to the path of the microwaves into the instrument and it is not influenced by pseudo- 
correlation. 

FEM bandpass measurements 

The most straightforward method for performing frequency response of a single FEM 
channel to the sky input is: 
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Figure 3.5.: Schematic of LFI radiometers highlighting pseudocorrelation strategy 



1. the FEM is placed in cryogenic environment 

2. a known stable input temperature source is attached to the sky load input 

3. a mixer and an sweeping frequency source is set up in order to selectively measure 
noise as a function of frequency 

4. stop the 4 KHz switching (otherwise the signal from a single channel switches every 
1 / 4096 s between sky and ref signal) 

5. set the phase switches polarisations in a configuration where the current channel 
receives the sky signal. Opposite Phase Switches configurations both have the 
same output correspondence, e.g. arm 1 phase switch direct and arm 2 ph/sw 
shifted by 180 is equivalent to 1-shifted 2-direct. In principle the Phase Switches 
have different attenuation in the abovementioned configurations, but this difference 
is indeed very small because they were balanced apply the right control currents 
before performing the tests. 

Therefore, Jodrell Bank Observatory, 30 and 44 GHz FEMs manufacturers, per- 
formed tests in only one of these configurations, while Elektrobit, manufacturer 
of 70 GHz RCAs, performed tests in both configurations. However, the difference 
between the two configurations is as average 0.2 dB, with peaks of 0.9 dB, the gains 



in both configurations are overplotted on the top of the RCA22 plot in figure 3.7 

For consistency with the data available for 30 and 44 GHz FEMs, for 70 GHz the 
mean of the two tests were taken as FEM response. 
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6. Measure the output signal with a Noise Figure Meter as a function of the signal 
frequency and amplitude 



7. Perform all the aforementioned steps at 2 different input temperatures; measuring 
the output at 2 different input power it is possible to compute the gain as ^W*^ 
and consequently the T no i se . 

For 70 GHz channels, the cryogenic tests consisted only in gain measurements, there- 
fore the noise temperatures used in the model were taken from the previous test cam- 
paigns on the LNA amplification stage, named ACA (Amplifier Chain Assembly), at 
ambient temperature. Integrated noise temperatures were measured in cryogenic con- 
ditions, therefore it is possible to normalise T no i se measured at ambient temperature in 
order to obtain the cryogenic ones when integrated. 

It is important to note that bandpasses, which are the main focus of QIMP, do not 
depend on either FEM or BEM T no i se because they are gain bandpasses, based on the 
proportionality between output and input, therefore not influenced by static contribu- 
tions: 



Vout = Gbem * (Gfem * {Tin + Tnoise FE M)) + VnoiseBEM (3.2) 
Vout = Gbem * G F em * T in + G B em * G FE m * Tnoise FE M + VnoiseBEM(3.3) 
V out = G*T in + Vnoise (3.4) 

It is clear that only FEM and BEM gains contribute to the G term of the last equation, 
which is the bandpass response, all noise terms converge to the Vnoise term. 

This modelling could be refined when the model will be extensively validated using 
the results of the LIS tests, which are tests specifically dedicated to understanding the 
integrated output volt signal behaviour as a function of input temperature. 

FEM at 30 GHz were measured just between 25 GHz and 35 GHz, therefore an 
extrapolation has been necessary in order to have a total simulated bandpasses frequency 
span between 21.3 and 40 GHz. Extrapolation slope is based on VNA measurements 
that were performed in warm conditions on a Qualification model FEM; its offset instead 



was based on the interface with the measurements. Figure 3.6 shows the result of the 
extrapolation on RCA 27. 

FEM Linear model 

The linear model fundamental FEM equation is: 

V out = G(T in ) + V sys (3.5) 

The output is the sum of the Gain multiplied by the input power and an offset which 
doesn't depend on the input temperature and which is due to the thermal noise emission 
from the electronic components. Usually instead of referring to the noise contribution 
as a volt offset, it is divided by the Gain and referred to as the System Temperature 
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Figure 3.6.: Example extrapolation at FEM level performed on a channel of RCA27 
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Figure 3.7.: FEM Gain and Noise Temperature as a function of frequency, an example 
from LFI27, LFI24 and LFI22, for LFI 22 gain and noise temperature are 
on different plots because they were measured independently 



(T sys ): it is the equivalent input temperature that would produce the same output in an 
ideal amplifier (i.e. a component without noise). 

V out = G(T m + T sys ) (3.6) 

It is very convenient because it can be directly compared to the input temperature and 
understand what is the relative importance of input signal and electronics noise in the 
output signal. During the data analysis Gain and System temperature are computed for 
each frequency step, typically 0.1 GHz, using the outputs at high and low temperature: 

1. Compute the Gain in dB as [20log(VoutT hi h — Voutx low ) / (ViriT high — Vinr low )] 

2. Compute the system temperature T sys extrapolating linearly the volt output until 
the intersection with the negative part of the volt axis. 

The next figures are examples of 30,44 and 70 GHz FEM gain and system tempera- 
tures. 
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Model 



Focusing on a single channel, and ideally freezing the 4 KHz switching when it receives 
the sky signal, it is possible to study statically the bandpass response. 

The purpose of QIMP is to characterise statically the radiometer properties: in this 
context, each single channel can be seen as an independent radiometer frozen in one of 
the phase switches configurations when it is looking to the sky signal. Actually, there 
are 2 opposite phase switches configurations where the same BEM channel is connected 
to the sky, i.e. the phase switches can shift the RF signal phase by 180°, changing the 
'sign' of the incoming radiation, if none or both the phase switches change the sign, the 
output correspondence to the input signal is the same. Due to the effect of hybrids, each 
channel bandpass response includes the contributions from the Low Noise Amplifiers of 
both FEM legs. 

In the QIMP schematic, the FEM has been modelled as an S-parameter component 
with gain given by the performance tests on each FEM channel and noise figure com- 
puted by the system temperature: F = WlogW(^fc + 1). FEM Sll and S22 were 
measured only at ambient temperature, in cryogenic conditions they are expected to be 
about -7 dB, therefore QIMP contains a value of -7 dB flat over the band. FEM Sll 
is not very important because OMT reflection is low, so the power reflected back negli- 
gible. S22 instead would be interesting because it influences the losses due to standing 
waves into the waveguides. For this reason in the future I plan to investigate if it is pos- 
sible to predict cryogenic S22 by properly scaling the coefficient measured at ambient 
temperature . 

Therefore the RCA is studied in an ideal configuration, which could never be obtained 
in the real hardware: all the 4 RCA channels are looking simultaneously to the sky 
signal, useful for characterising at once the broadband behaviour of the instrument. The 
main disadvantage of this modelling approach is that it is not possible to study isolation 
between sky and reference signals, which needs the full pseudocorrelation design to be 
implemented. 

Notice that the response to the reference load input is in principle different from the 
response to the sky, because the receiver is asymmetric upstream the first hybrid. In 
fact the sky channel is connected to the feed horn - OMT assembly and the reference 
channel is connected to the 4K reference horn which is connected directly to the 70 GHz 
FEMs or by means of 2 short waveguides to the 30 and 44 GHz FEMs. 

The effect of this asymmetry on RCA25 was analysed by Andrea Catalano using 



QIMP, the results, in Figure 3.8 show that the difference between sky and reference 
bandpasses is below 0.15%. I extended this analysis to all LFI channels confirming that 
the impact of the receiver asymmetry on bandpasses is negligible. 

3.3.6. Waveguides 
Description 

Waveguides link the 20 K Front End Module to the 300 K Back End Module, and are 
mechanically and thermally connected to the 3 V-Grooves at 60, 80 and 140 K which 
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Figure 3.8.: Comparison between sky and reference bandpasses on RCA25 
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thermally decouple the warm stage from the cold. 

Waveguides electromagnetic behaviour is strongly influenced by the large temperature 
gradient. 

Each waveguide is made by a twisted copper section connecting the FEM to a straight 
stainless steel that thermally decouple the 20 K from the 50 K stage. Then a gold plated 
stainless steel straight section connects this stage to the 300 K BEM. 

Frequency response measurements 

IFP (Istituto di Fisica del Plasma) tested all the waveguides at ambient temperature 
characterising Return and Insertion losses with a Vector Network Analyser, details about 
the measurements are available in |10j . 

Because measurements were performed with impedance-matched connections at the 
waveguides ends, they are not representative of the standing waves arising when the 
waveguide is connected to FEM and BEM and therefore are source of important impedance 
mismatch. 

For this reason QIMP implements an analytical waveguide simulator built using many 
sections of a rectangular waveguide component and ideal transformers. The waveguide 
model was validated against measurements performing simulations at ambient temper- 
ature, after validation the model temperatures were modified in order to simulate the 
waveguides in nominal conditions, i.e. with the nominal temperatures at FEM, V-grooves 
and BEM connections and linear interpolation between them. 

Rectangular waveguide device implementation in QUCS 

In summer 2008, under my request, Bastien Roucaries, one of QUCS developers, im- 
plemented a rectangular waveguide component model, which was not present before in 
QUCS. 

The model is based on the analytical waveguide treatment in appendix [B] 
The rectangular waveguide model was implemented using a S-parameters component. 
Its implementation has been quite straightforward thanks to the modular conception of 
QUCS. The code is constituted by about 600 lines including extensive comments and 
GUI part. The core part was done reimplementing six virtual functions of the component 
base class. 

Upon the model developed by Bastien Roucaries I added the embedded computation 
of the resistivity of gold, stainless steel and aluminium based on device temperature. 
This feature was implemented using different empirical formulae in different temperature 
ranges. 

The waveguide component was validated against Agilent ADS results, using a reference 
waveguide model as follows: 

• single rectangular waveguide component 

• 5.69x2.8 mm size (same as 44 GHz LFI waveguides) 

• 800 mm length 
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Figure 3.9.: Comparison of S21 [dB] of a gold plated rectangular waveguide as a function 
of frequency 



• material is gold 

• simulation temperature is 20 °C 

Both S parameters and noise simulations voltage results were compared and the match- 



ing is very good, Figure 3.9 shows S21 as a function of input frequency of ADS and QUCS 
overplotted. 



Model 

Waveguides were modelled using many sections of the rectangular waveguide component, 
which consists in an analytical model based on the propagation of the Transverse Electric 
mode 1,0 (TE10). This was first implemented in ADS and then ported to QUCS. 

The transmission line has a reference impedance of 50O, while waveguides have a much 
higher impedance. Interfacing the waveguides directly to the transmission line would 
lead to strong losses, which are not present in the real hardware, because all components 
are physically very well matched. 

Therefore ideal transformers were used for the purpose of simulating the good match- 
ing between waveguides and the rest of the RCA, and only a small mismatch of -35 dB 
was left, as found during hardware measurements. 

The waveguide model is composed by 2 sections each enclosed by 2 ideal transformers: 



copper section, see figure 3.10 



stainless steel + 3 sections of gold plated stainless steel (figure 3.11) modelled by 



groups of single waveguide components at different temperatures. In figure 3.12| we 
show the detail of the stainless steel waveguides modelled joining many waveguide 
elements each at a different temperature interpolated linearly. 
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Figure 3.10.: Snapshot of the QUCS implementation of the copper waveguide section 
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Figure 3.11.: QUCS schematic of the group of stainless steel and gold plated sections 
enclosed by the ideal transformers 
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Figure 3.12.: QUCS schematic of the implementation of the stainless steel section, 
using many rectangular waveguide components at a linearly changing 
temperature 



Thanks to this approach it is possible to take into account the influence of temperature 
gradient on physical temperature and on material resistivity (for more details please refer 
to Paola Battaglia's work [3]). 

This waveguide model was tested by Cristian Franceschet against the Flight hardware 
tests using ambient temperature, the results showed a very good matching between 
measurements and simulations [16 1. 

As already mentioned, simulations were performed at ambient temperatures only for 
validation against the measurements; for producing the complete RCA bandpasses, the 
nominal temperatures are used, which are ~ 300-ftT at the Backend, ~ 140K, ~ 80K and 
~ 60K at V-grooves interfaces and ~ 20K for the entire copper waveguide. 

3.3.7. Back End Module 
Description 

The Back End Module receives a Microwave signal from the waveguides and outputs 
a DC volt signal that is then integrated in time and digitised by the Data Acquisition 
Electronics (DAE). Each of the 4 channels of a Back End Module contains: 

• a Low Noise Amplifier RF amplification stage 

• a bandpass filter 
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Figure 3.13.: BEM input power to output tension transfer function in dB at high (48°C), 
nominal(26°C) and low5°C temperature 

• a detector diode which integrates the RF signal to a DC output 

• a DC amplifier 

Frequency response measurements 

The objective of the bandpass characterisation of the BEM was to correlate the DC 
output (V) to the power RF input (W) as a function of the input signal frequency. The 
DC output in Volts comes from a diode and therefore it is proportional, considering a 
BEM linear behaviour, to the input power plus a noise term. 

The tests were performed by introducing a monochromatic input wave of known power 
(for example for 30 and 44 GHz was -60 dBm, which is the typical power entering the 
BEM when observing Tsky ~ 3K) sweeping over the interesting bandwidth. 

BEM response can be modelled as: 

V ut = GbEM * Wi n + Vnoise (3.7) 
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Out of band the Gain is about zero, therefore the output is a good estimation of the 
BEM noise voltage; therefore the BEM gain can be measured by setting the out of band 
output voltage as zero level and computing the transfer function between output voltage 



and input power, see figure |3.13| for LFI27. 

The logarithmic V over W gain without the offset and the BEM noise temperature 
(which works as an offset independent from the frequency) characterise the BEM fre- 
quency response. 

Model 

The RF signal coming from the waveguides is summed with a noise generator based on 
the measured BEM noise temperature and then amplified frequency by frequency by the 
V/W transfer function. The output is a bandpass response of the complete RCA which 
can be then integrated over the frequency in order to have the output in Volt that can 
be compared to the hardware output volt. 

BEM non-linear model During the Flight hardware test campaign, 30 and 44 GHz 
RCAs showed signal compression at BEM level, i.e. the BEM volt output doesn't increase 
linearly with the input temperature. This effect has a big impact on noise temperature 
computation, because noise temperature is computed by linear fitting volt outputs at 
different temperatures, computing the gain and the offset. 

The intercept with the negative part of the temperature axis is equal to the noise 
temperature; if the radiometer response is not linear, but we are assuming that it is 
linear, its slope is underestimated and its noise temperature is overestimated. 

Therefore, Fabrizio Villa and Luca Terenzi (|35j) applied a non linear analytic model 
to LFI radiometers that was used to fit the test data in order to evaluate the overall 
RCA Gain, System temperature and the BEM compression factor. 

The model links the input temperature, Ta, to the output voltage V ou t through the 
following relationship: 

Vout = Gun l + b.G lm \T A + T N) ) • {TA + Tn) (3 - 8) 

where Gu n is the uncompressed gain, b is the compression factor, Ta is antenna tem- 
perature and T/v noise temperature. The compression factor is a single scalar value per 
channel which represents the BEM Gain dependence on the input power. The non linear 
behaviour arises from a compression effect both in the RF amplification stage and in the 
diode. 

Considering that the BEM response is non linear, it is necessary to admit that the 
measured data we have from the performance campaign are already affected by non 
linearity and we need first of all to compute an equivalent linear transfer function, Gu n 
from these data. 



Using equation 3.9 it is possible to compute the linear gain Gu n knowing the test 
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input power, the measured compression factor b. 



G lm = ^ (3.9) 

The BEM non- linear model is not yet been implemented in QUCS, but it could be 
easily performed. However, due to the fact we do not have extensive measurements 
of BEM frequency response with different input power level, it is still under discussion 
whether the effect of the compression on the bandpass is a flattening, a negative offset or 
a combination of both. Therefore it is necessary to implement both models into QIMP 
and deeply study their impact on bandpass response. 

QUCS implementation Due to a bug in the ADS software related to a cascade of 
multiple noisy components, the Advanced LFI RF model implementation of the BEM 
was very complex and more difficult to manage and improve. 

In QIMP instead the implementation is by a single S-Parameters device; its return 
loss as a function of frequency and its integrated noise temperature must be derived 
from BEM measurements. 

Because measurements provide output voltages versus input power, so that: 

Gbem = 777^ (3.10) 

we need to manipulate the previous equation in order to convert the power measure- 
ments to amplitudes, i.e.: 

Gbem = ^ (3.11) 



V Gbem — ^Z^L = (3.12) 



VVo^t = VGbem • V50V in (3.13) 

where the factor 50 is due to the standard 50f2 transmission line impedance. 

Therefore we have a simple linear relation between the input voltages and the square 
root of the output voltage. 

It is therefore just necessary to create an equation object and apply the square to have 
the output voltage Vout- 

3.3.8. QIMP simulations 

The QUCS simulation is run on the complete schematic defined following the previous 



sections, see figure 3.14 



The result is a voltage output in Volts as a function of frequency; the simulations 
approach is very similar to swept source measurements where a monochromatic input 
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Figure 3.14.: Complete schematic of the LFI Radiometer model implemented in QUCS 



power is swept through the bandpass and the BEM voltage output is recorded; the 
main difference is that in the test there is also a small contribution to the input due 
to the background temperature of the sky load all over the receiver band. However, 
this contribution is negligible because the injected microwave signal is about 3 orders of 
magnitude higher than the background temperature. 

All the RCA bands are swept with a 0.1 GHz step: for each step the power input to 
the Feed Horn is a rectangular step (in the frequency domain) centred in the sweeping 
frequency and 0.1 GHz wide, with the power level determined by the input black-body 
thermal emission spectrum. 

Taking into account the frequency response of each component and multiple reflections, 
including the waveguides model, the LFI Radiometer model computes how the input 
frequency spectrum is modified by the effect of each component. 

An example is provided in figure 3.15 that shows plots of the bandpass along a 30 
GHz RCA: it is important to note that neither the bandshape nor the total power in 
any point along the RCA before the BEMs have been measured experimentally during 
the test campaigns. 

The first objective of the model port to QUCS was to reproduce exactly the same 
results concerning output bandpasses; after the waveguide model was implemented into 
QIMP, the result was reached successfully, in the figures 3.16 and 3.17 I show an exam- 
ple of the comparison between ADS model, QUCS model and measurements for each 
frequency 
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Figure 3.15.: Bandpasses shown along the RCA 
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RCA20 channel RIDO 




Figure 3.17.: Comparison of ADS, QUCS and measurements bandpasses, LFI20S-10 



3.3.9. Automation 

QIMP includes also a set of Python programs needed for automating most of the pro- 
cedure of running simulations. I will just review their main features. 



They are based on the software package SciPy, see http : //www . scipy . org a general 
purpose scientific and numerical software for Python, in particular on: 

• NumPy: multidimensional array package 

• Matplotlib: 2D plotting library 

• Ipython: interactive Python console 

The programs are implemented following Object Oriented Programming and are stricly 
independent, each stage runs separately by reading inputs and writing results on stan- 
dard ASCII files, the stages are: 

ADS2QUCS: Format conversion from ADS to Touchstone, merging of new input data, quick 
modification of input data 

lfisim: Batch simulation run using qucsator, QUCS GUI is never launched, and export 
of bandshapes from QUCS to ASCII format; a simulation of all 44 channels takes 
about 5 minutes. 

ba_lib: Bandshape analysis tool, for batch processing and interactive analysis: 

— gain bandshape from simulations at 2 different temperatures 

— normalization 

— plotting 

— comparison with Swept source measurements and ADS simulation 

— export to text format 

Moreover, thanks to Ipython, it is easy to switch from batch processing to an interac- 
tive data analysis session which is very useful when investigating a specific issue. 
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3.4. Frequency response measurements 



The most important validation of the QIMP model relies on the comparison between 
the modelled and measured bandpasses. 

This section contains a description of the measurements setup and results which are 
important for understanding the results of the comparison. 



3.4.1. 30 and 44 GHz Radiometers 
Experimental setup 

30 and 44 GHz RCAs were assembled and tested in Thales Alenia Space in Vimodrone 
during Spring 2006; each RCA was tested separately at ambient and cryogenic temper- 
ature. By requirements, the aim of the test was to measure relative frequency response 
of the instrument, i.e. what is the response relative to its own maximum. 

The bandpass was measured by injecting a calibrated monochromatic source into the 
sky horn sweeping through the operational band and recording the DC output as a 
function of the input frequency, see figure |3.18[ 
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Figure 3.18.: Schematic of the experimental setup for 30/44 GHz Swept Source Tests: a 
monochromatic input microwave is injected into the sky load, its frequency 
changes with time through the radiometer band. From the recorded output 
voltage the bandpass has been determined using Eq. 2.16 
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The Agilent synthesised microwave generator HP83650B was located at ambient tem- 
perature outside of the cryogenic chamber; its output on RF coaxial cable was transferred 
to rectangular waveguide (WR28 for 30 GHz and WR22 for 44 GHz) by an Agilent tran- 
sition. The transition was connected with a short waveguide section to the cryogenic 
chamber by a RF vacuum feed-through with a kapton window. Inside the chamber, 
another straight stainless steel waveguide section acts as a thermal break, to reduce 
thermal inflow from the outside environment; the last section of the transmission line is 
a 1-meter long flexible waveguide which allows the connection to the sky load. 

The sky load is a cylindrical black cavity 20cm high and with a diameter of 20cm; its 
back wall, positioned in front of the horn, is a bed of Eccosorb pyramids 5 mm large 
and 30 mm high. Sky load is composed by modular absorbing panels screwed to an 
aluminium plate, they are free to float on it in order to compensate for different thermal 
contractions. 




Figure 3.19.: Picture of the injector connected to the sky load for a 30 GHz RCA; the 
large hole in the centre of the sky load cylindrical base is dedicated to the 
feed horn 



The injector, see figure 3.19 is a open ended waveguide located on the same side 
of the horn with the axis parallel to horn axis; therefore the signal was reflected into 
the horn by the sky load back wall. The sky load return loss is about 60 dB, so that 
with the input power set to -35 dBm, the expected power delivered to the Feed Horn 
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is about -95 dBm. Not knowing exactly the input power permits only to compute a 
relative bandpass, which is coherent with test requirements, but forces the comparison 
with modelling to be based just on normalised bandpasses. 



Results 

The measurement campaign gave poor results for 2 channels of the same polarisation 



axis of each 30 GHz RCA due the presence of 1 GHz spaced ripples, see upper figure [3^20 
and a very low and noisy Volt output. Investigation at warm temperature showed that 
the power delivered to the Feed Horn was highly dependent on injector orientation; 
the ripples disappeared by connecting the injector with its polarisation plane at 45 
degrees with respect to OMT Main arm and Side arm polarisation planes (refer to [36 ). 
Unfortunately it was not possible to repeat the tests in cryogenic chamber, therefore 
frequency characterisation of these channels can only rely on simulations. 

Among the eight 30 GHz channels, four of them show ripples as highlighted in the 



previous paragraph, while the remaining 4 are very clean, see lower figure 3.20, and 
extends down to 26.5 GHz; this limit was established considering the nominal bandwidth 
of 20% of the central frequency, e.g. 6 GHz [27-33 GHz] for 30 GHz and that 26.5-40 
GHz is the nominal band of WR28 waveguides, used for injecting the signal into the sky 
load. However, the response at 26.5 GHz is still high and useful bandwidth is expected 
below this limit, but the shape of this low frequency cut was not measured. 

All twelve 44 GHz channels measurements are very clean and frequency covered range 



is sufficient, see figure 3.21 the response at 50 GHz is still higher than background noise 



but already one order of magnitude less than central band gain. 
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Figure 3.20.: Swept Source bandpass measurements of 30 GHz: RCA LFI28S-11 is an ex- 
ample of channels showing ripples due to injector orientation, while LFI27S- 
11 an example of a clean measurement, bandpasses measured are just rel- 
ative, therefore the normalisation is arbitrary, in this case the integral was 
set equal to the corresponding modelled bandpasses - see section 3.5 
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Figure 3.21.: Swept Source bandpass measurements of LFI25S-10, 44 GHz, bandpasses 
measured are just relative, therefore the normalisation is arbitrary, in this 
case the int egral was set equal to the corresponding modelled bandpasses 
- see section 3.5 
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3.4.2. 70 GHz Radiometers 



Experimental setup 

70 GHz radiometers were produced, assembled and tested at Elektrobit Ylinen in Fin- 
land, see [20] . The Swept Source test setup was different than for 30 and 44 GHz tests: 
the Microwave Generator was connected to a small sky load directly in front of the horn, 



see figure 3.22 




Figure 3.22.: Picture of the two 70 GHz horns connected with the small sky loads pro- 
vided by Ylinen, input signal waveguides and injectors are located on the 
left 



Thanks to this design it is possible to estimate the power effectively delivered to the 
horn by using a calibrated input source and applying corrections due to the connection 
waveguides and the injector; in this conditions the skyload delivers all the input power 
to the feed horn. 

The main issue of this system is the return loss of the sky load; the requirement was 



-30 dB, but the measured performance didn't meet the requirement, see figure 3.23 
The favoured hypotheses is that the observed ripples in the bandpass measurement are 
caused by strong standing waves between the feed horn and the sky load due to its high 
return loss. 



Results 

The standing waves had a big impact on measured bandpasses, which suffered of regularly 
spaced ripples between 1 and 2 GHz wide and with a peak to peak amplitude between 
5 and 15 dB. 
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Figure 3.23.: Measured return loss of the sky load used for testing the 70 GHz RCAs, 
the plot shows two different type of sky loads, #2, which was used for the 
tests, at different temperatures. 



The ripples are very similar in all 70 GHz channels, see an example in figure 3.24 



moreover ripples pattern is not regular enough for removing it with software manipula- 
tion of the data without the risk of adding additional errors. 
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Figure 3.24.: Swept Source bandpass measurements of LFI19S-10, one of the 70 GHz 
channels, bandpasses measured are absolute 
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3.5. Comparison of measurements and simulations 



As written in section 3.4, 30 and 44 GHz RCA Swept Source measurements provided 
only relative bandpasses while 70 GHz tests provided absolute bandpasses between input 
power and output voltage. 

Therefore, for 70 GHz the comparison process is straightforward, simulations bandpass 
is computed by: 

q _ Phigh ~ Plow ^\ 
(Vhigh — Vlow) 

where the subscripts refer to simulations performed at 2 different temperatures. The 
dimensions of the result is [W] / [V] and can be directly compared to Swept Source Tests 
results as explained in section |3.4.2| 

For 30 and 44 GHz channels the input power is very high compared to radiometer 
noise temperature, therefore T no i se effect on the output is negligible and it is possible to 
compare the normalised bandpass directly with the normalised simulated gain. 

In order to compare measurements to simulations measured data have been normalised 
as model results: 

(3.15) 

where: (3.16) 
A = f™ as{V ] • ( G model du (3.17) 



J G meas du 

Comparisons of measurements and simulations for all LFI channels are plotted in 



figure 3.35 and figure 3.36 at the end of the chapter. 



3.5.1. 30 GHz Radiometers 

One of the noticeable results of the QIMP simulations has been to provide information 
about the bandwidth below 26.5 GHz, down to 21.3 GHz. In fact RCA measurements 
were available only from 26.5 GHz to 40 GHz. 

FEM response at lower frequencies has been obtained extrapolating the available data 



using as a reference VNA measurements performed on a larger bandwidth, see 3.3.5 



OMT response instead was extrapolated by exploiting its similarity with 44 GHz OMTs, 



which were measured on a much larger frequency span, see 3.3.4 



At the end of the chapter, in figure 3.36, we show the comparison between simulations 
and measurements for all the 30 GHz channels. 

It is important to note that measurements convey in this case, and for 44 GHz as well, 
only a relative information. They have a free normalisation factor which translates into 
an offset in logarithmic scale; therefore they were normalised using the model output as 
a reference. 

In general there is good qualitative matching between measurements and simulations, 
which supports the hypotheses that there are no major influence among the components, 
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Figure 3.25.: 30 GHz channel LFI28S-11, an example of the polarisation issues on half 
of the 30 GHz channels: measurements were normalised considering model 
output as a reference, bandpasses units are dB : 10 log 10 ([V]/ [W]) 
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Figure 3.26.: 30 GHz channel LFI27S-11, good agreement between measurements and 
simulations 
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Figure 3.27.: 30 GHz channel LFI28M-01, general good agreement beside a discrepancy 
in the central part of the band 
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i.e. the performance of radiometer devices doesn't change significantly when assembled 
into an RCA. 

On the other side this is also a good indirect test of performance tests at component 
level, which have been effective in characterising radiometer frequency response. 

Some channels show strong discrepancies due to issues on the measurements setup 
while others show smaller discrepancies still under investigation. 



In particular Figure 3.25 shows an example of a channel affected by high ripples due 



to the measurements test setup, see section 3.4.1 in this case comparison is useful just 



for fitting the high frequency gain drop slope. Channels affected by this issue are 4 on 
the total of 8 channels. 



Figures 3.26 and 3.27 shows comparison from RCA 27 and 28, these channels are very 
similar to their own twin channels; discrepancy above 36 GHz is due to the thermal noise 
of the BEMs, which is present in the measurements but not taken into account by the 
model. 

RCA 27 comparison is very good while RCA 28 shows an important discrepancy in 
the central part of the band which is still under investigation. 



3.5.2. 70 GHz Radiometers 

Due to the different test setup, 70 GHz measurements gave absolute bandpasses, there- 
fore the comparison with the model is both on the output level and on the band shape. 

In this case measurements are slightly broader than simulations, 57.5-82.5 GHz versus 
60-80 GHz, but they are highly affected by ripples very likely caused by standing waves. 



As highlighted in section 3.4.2, the main issue concerning 70 GHz Swept Source Tests 
is the strong presence of ripples in the bandpass; the modulation of the superimposed 
ripples pattern is even larger than the modulation of the radiometer bandpass, its re- 
moval, indeed, is not straightforward and may introduce uncontrolled systematics in the 
result. 

However, the difference between model and measurement looks really consistent with 
a standing wave pattern; a good evidence supporting this hypothesis is that Elektrobit 
performed a preliminary test of FEM and BEM alone with direct signal injection in 



waveguide and its result doesn't show traces of ripples, see figure 3.28 



Most of the channels show high consistency between model and measurements (see, 
e.g., top of figure 3.29); few channels have an offset in logarithmic scale, which means 



they differ by a multiplicative factor, see middle of figure 3.29 The cause of this effect is 



not clearly understood, but its impact on bandpasses modelling is low because the main 
objective of bandpass modelling and simulation is just to deliver relative bandpasses. 
Absolute level consistency is good for increasing confidence on radiometers modelling 
but it is not a required feature. 



RCA 23 instead, see bottom of figure 3.29 is a peculiar case, because all its 4 channels, 



in the centre of the band, show a model bandpass higher than measurements bandpass; 
this issue is under investigation, the favoured hypotheses is a problem on model input 
data. 
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Figure 3.28.: Ylinen Tests on FEM and BEM with signal injection in waveguide directly 
into FEM 
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Figure 3.29.: 70 GHz channels LFI21S-10 on top, LFI22M-00 on bottom, bandpasses 
measured and simulated, for 70 GHz channels the comparison is on absolute 
value and not only on normalised shapes 




Figure 3.30.: 70 GHz channel LFI23M-00, bandpasses measured and simulated, for 70 
GHz channels the comparison is on absolute value and not only on nor- 
malised shapes 
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Figure 3.31.: RCA 25 44 GHz channel S-ll bandpass measured and simulated: measure- 
ments were normalised in order to have the same integral as the model, a 
systematic discrepancy between 44 and 49 GHz is present in all 44 GHz 
channels and it is still not understood 



3.5.3. 44 GHz Radiometers 



Figure 3.31 shows an example of a 44 GHz channel; few notes about the comparison: 

• measurements show no systematic effects 

• model and measurements have the same frequency coverage 

• measurements were normalised to the same integral of simulations, using the same 
procedure applied on 30 GHz channels 

• agreement is good apart from a systematic discrepancy in all 44 GHz channels 
from 44 to 49 GHz, where simulations have an higher response 

The discrepancy is important and at the moment it is not clear what is the cause of 
this effect, among the hypotheses there are both model deficiencies: 
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• systematic error in a subcomponent campaign test 

• insufficient modelling of the FEM return loss on the back side (facing the waveg- 
uides) 

or measurements issues: 

• sweeper input power loss 

• error on the measurements data analysis process 

The last possibility at present is the favourite because a preliminary reanalysis on the 
data acquired during swept source tests on just one channel gave a result much more 
similar to QIMP simulations that the first analysis. However, this is not true for all 
channels, therefore I plan to analyse all the data using a different software in order to 
investigate this issue further. 



3.5.4. 30 GHz Spare radiometer 

In 2008 Stefania De Nardo, during her master thesis [12], performed new bandpass mea- 
surements in warm conditions on a 30 GHz spare radiometer. Spare units are equivalent 
to flight units and were built to be used for replacing flight units in case of damage. 

These tests were performed on a much larger band, from 21.3 to 40 GHz, and were 
also performed both on the complete RCA and on different chain sections. In particular 
here is the list of tests performed, in which the signal was injected at different points: 

• into the skyload, with 3 different input waveguide orientation 

• into the OMT, connecting the injection waveguide directly to OMT main and side 
channels 

• into the Waveguides 

• into the BEM alone 

In all cases the DC signal output from the BEM was sampled. 

Tests with skyload injection suffer strong ripples probably due to the presence of 
standing waves in the skyload, tests with injection into the OMT, instead, are very 
clean. The feed horn effect on bandpass is very low, because its losses in general are 
negligible, therefore tests with injection at OMT level can be regarded as representative 
of the RCA response. 



Figures 3.32 and 3.33 show the comparisons between measurements and simulations 



for main and side channels: 

• the agreement is very good, even below 25 GHz, where FEM and OMT data were 
extrapolated in the model, and this is an important validation of QIMP results 
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Figure 3.32.: Comparison of measurements and QIMP simulations on 30 GHz spare 
MAIN channels 



• in Main channels the discrepancy is mainly in the dips at the centre of the band; 
this could be due to the power reflected back by the FEM which is reflected by 
the test setup again into the FEM. Such an effect tends to flatten the band and 
in nominal configuration it is not present because the skyload absorbs the power 
reflected back by the FEM. 

• in Side channels the effect of bandpass flattening is weaker, probably because it is 
already flatter than Main band, however there is a discrepancy between 33 and 35 
GHz which is still under investigation. 



3.6. Conclusion 



Swept Source Measurements suffer many systematic effects, mainly standing waves on 



70 GHz (section 3.4.2) and polarisation issues plus low frequency coverage on 30 GHz 



(section 3.4.1). These effects are zero order effects and are so important that prevent 
using SPR as best bandpass estimator, and make the mathematical difference between 
measurements and simulations unmeaning for judging their likelihood. 

On the other hand single components performance tests (which are less prone to 
systematic effects because of simpler setups) and QIMP simulations provide the best 
knowledge of LFI radiometers frequency response. 
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Figure 3.33.: Comparison of measurements and QIMP simulations on 30 GHz spare SIDE 
channels 



The QIMP model has been successfully validated against the bandpass measurements. 
Some discrepancies remain at 44 GHz, this case is still under investigation. 

The best estimation of Planck LFI radiometers bandshapes up to now is represented 
by the QIMP simulations, obtained by exploiting single components performance tests. 
Simulations have still margin for improvement, thanks to the analysis of spare com- 
ponents tests. For the first time they included SPR tests at ambient temperature of 
sections of RCA in order to better understand how the microwave signal is transmitted 
by each component. 



3.6.1. Use of bandpasses into data analysis 
Temperature measurements 

CMB anisotropics signal is highly contaminated by foregrounds emission and data anal- 
ysis strategies rely on their different frequency dependence in order to distinguish them 
from the CMB. 

Figure 1.2.4 in chapter [T] shows the different frequency dependence of dust, syn- 



chrotron, free-free emissions and CMB; thank to its wide frequency coverage, Planck 
will be able to separate the different components by combining maps from different 
channels. 

Both radiometers and bolometers, however, have a broadband response, therefore their 
output is a convolution of the input spectra with their own frequency response. 
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A fine correction can be applied at map level by an iterative process which separates the 
different components thanks to their spectral indexes, convolve them one by one with 
the instrument frequency response, and compares their sum with the detected map. 



Polarization measurements 

The new frontier of precision cosmology is CMB polarization measurements; several HFI 
channels at 100, 143, 217 and 353 GHz are sensitive to polarization thanks to Polarization 
Sensitive Bolometers, |22j . 

Each LFI Feed Horn, instead, is connected to an OrthoMode Transducer which splits 
the orthogonal components to be detected by two twin radiometers. Each radiometer is 
calibrated separately against the CMB dipole signal and polarization measurements are 
performed differencing the signals of the 2 radiometers connected to the same Horn. 

The main issue is that foregrounds have different spectral indexes with respect to the 
CMB dipole which determined the calibration, therefore spurious polarization signal is 
generated when differencing channels with a different bandpass, producing leakage of 
foregrounds temperature signal to polarization. 

Using the radiometers frequency response, it is possible to estimate the impact of 
bandpass mismatch on CMB maps and power spectra, this work was performed for 70 
GHz radiometers by Torsti Poutanen in [29j. 



3.6.2. Delivered bandpasses 

QIMP bandpasses convey information both on relative bandshape and absolute gain; 
however, there is a more reliable technique for measuring the absolute gain, which con- 
sists on dedicated tests which measure the Volt output as a function of input power and 
then fit radiometer response with an suitable model, details in |15j . 

Therefore bandpasses need just to characterise a relative response, and they are nor- 
malised: in detail a weight W{ is computed by integrating the bandpass G(v) over 0.1 
GHz, they are then normalised to have a sum of 1. 

Absolute frequency response can then be computed by multiplying QIMP bandshapes 
by the absolute gain, called photometric constant. 

Tabulated bandpass files are available at LFI Data Processing Centre in Trieste or 



from Andrea Zonca at request upon approval of the LFI Instrument Team. Figure 3.34 



shows all the 44 channels normalised bandpasses on the same plot in logarithmic format. 
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Figure 3.34.: All the 44 LFI channels bandpasses released on September, the 18th, 2007, 
bandpasses are normalised to unit integral and plotted in logarithmic for- 
mat computed as 10 log Wi where each Wi is the radiometer gain G{v) inte- 
grated over 0.1 GHz, the separation in frequency, and normalised to sum 
equal to 1. 
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3.7. Comparison of all the channels bandpasses 

The following figures plot channel by channel all the comparisons between the measured 
and QIMP bandpasses. 
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Figure 3.35.: 70 GHz RCAs comparison between measurements (black) and model (blue) 
results 
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4. Thermal aspects 



4.1. Introduction 

Cryogenics and thermal issues were a strong driver in the design of the Planck Satellite, 
HFI bolometers need a temperature of 0.1 K to work, while LFI radiometers need a 
Front End temperature of ~ 20K to achieve the required sensitivity. 
Planck cooling system is composed by, see figure |4~T| 

• three passive radiators at 140, 80 and 60 K which decrease gradually the temper- 
ature from the 300 K of the Service Module, to the interface with the first cooling 
stage 

• the Sorption cooler is the most important cooling system onboard Planck because 
it cools the complete focal plane, including both LFI and HFI, from 60 to about 20 
K. Technically it is a vibrationless hydrogen cooler through which the hydrogen is 
pumped exploiting a chemical sorption process. It is composed by six compressors, 
each compressor is periodically cycled between heating and cooling phases, each 
phase lasts 667 seconds and a complete cycle 667s x 6 = 4002s. It is capable of 
providing ~ 1W of cooling power at 20 K to cool the LFI radiometers and as 
precooling stage of the 4K. 

• the Joule-Thomson cooler is a mechanical helium compressor, it cools the HFI 
external shell to 4K. The 4K shell is very important because the it includes both 
the HFI focal plane, where its horns are located, and the LFI reference loads. In 
order to improve the stability of the 4K stage top plate, where horns are located, 
a PID controller commands a heater. LFI 70 GHz reference loads are very near 
to the top plate and they benefit from the PID action, 30 and 44 GHz are located 
further and their stability is lower. 

• the Diluition cooler is used by HFI to keep the bolometric detectors at 0.1 K 

4.2. 4K box thermal transfer functions 
4.2.1. Thermal transfer functions 

Before discussing the work performed on data analysis of the measurements, I'd like to 
introduce the concept of thermal transfer functions. 

If a solid body has a boundary imposed temperature which is sinusoidally fluctuating: 

T(x , t) = T + T s sin{2-Kut) (4.1) 
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Figure 4.1.: Schematic of Planck cooling chain 

where To is the background temperature, T s the amplitude of the sinusoidal fluctuation 
and v its frequency, the generic solution at any other point of the solid body is still a 
sinusoid: 

T(x, t) = T + T s 7(x, v)sin(2nvt + </>(x, v) (4.2) 

where 7(x, v) is the amplitude transfer function and (f>(x, v) is the phase transfer 
function. This derives from the fact that the equation of heat propagation is linear. 

The knowledge of the transfer functions for each frequency v allows to compute the 
fluctuations at a point x due to any boundary condition by working in the Fourier do- 
main: any smooth source fluctuation can expanded as a series of sinusoids, each sinusoid 
is transformed to the output by adding the phase and multiplying by the amplitude 
transfer function and the fluctuation at the new point is reconstructed by merging again 
the transformed sinudoids. 

4.2.2. Introduction 

This section presents the analysis performed on the data acquired during the thermal 
transfer functions measurements, performed at the Liege Space Centre in cryogenic 
chamber on the Planck Cryogenic Model in 2005. The Planck Cryogenic Model con- 
sisted in: 



the qualification model of the satellite, built for qualification and testing purposes 
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• HFI qualification model (CQM), the complete HFI instrument built for testing and 
qualification, just some redundancies were missing with respect to the flight model 

• LFI instrument was not involved in the test, it was replaced by a dummy structure 
needed for simulating its passive thermal behaviour 

During the test periodic thermal fluctuations with periods of 100, 666 and 4000 seconds 
were forced at LVHX1 and the resulting effects measured with the HFI thermometers 
on lower part of the 4K box (Cernox), on the side, where LFI 30 and 44 GHz reference 
loads are located (Ther4KLl,2) and on the top 4K bolometers plate (PID4C,R,N). 

The tests for 100 and 666 seconds were performed both with PIDs controllers on and off 
while the 4000 seconds only in PID on condition, due to the more important drift effect. 
In order to estimate the thermal transfer functions, an analysis has been performed in 
the time domain by computing the ratio between the fluctuation amplitudes at the input 
and at the output. 

This method is straightforward but it is not very reliable when the input and output 
signals are not perfect sine waves, therefore it was applied, for validation purposes, only 
in analysing data of the 666 PID off and 100 PID off in order to compare its results with 
an analysis in the frequency domain. 

Analysis using the Fourier Transform 

The thermal transfer function versus input frequency was evaluated by applying the 
Fast Fourier Transform both to the input and to the output time streams, after having 
removed the transient. The results were divided by their own sampling frequency in 
order to compare different type of sensors (for example LVHX1 has 1 Hz sampling while 
PID4* thermometers have 160 Hz). 

Finally the amplitude transfer function was computed by taking the ratio between the 
peaks of the output on the input FFTs, at the stimulus frequency; while the phase was 
computed making the difference between the phase of the output and the input peaks. 

4.2.3. Test setup 
Description 

The HFI CQM was integrated on the Planck CQM satellite. The Sorption cooler Pipe 
Assembly Cold End (PACE), the 4K cooler CQM and the dilution cooler CQM were 
also integrated on the satellite and connected to the instrument. 

The whole satellite was integrated in a vacuum chamber cooled at 20 Kelvins. 

The PACE was fed by Hydrogen bottles instead of compressor. 

The thermal transfer function measurement principle for the HFI at CSL was to 
generate various thermal stimuli at the 18K cold ends of the PACE by modulating the 
hydrogen flow in order to quantify the impact of in-flight Sorption Cooler fluctuations 
on HFI stages stability. 

Various pressure changes sequences in the PACE allowed to generate different periodic 
fluctuations of LVHX1 and LVHX2 with various amplitudes. From the result of the 
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various PACE fluctuation tests it is then possible to deduce the thermal transfer function 
between LVHX1 and the other HFI stages. 

The complete test report is available from |18j . 



Thermometers location and performance 



Figure 4.2 shows the location of the thermometers on the HFI 4K shell and includes a 
table with the sampling frequency of all the thermometers. Thermometers are divided 
in three groups: 

• high precision and high frequency thermometers share the same acquisition elec- 
tronics of the bolometers, due to their key role for the processing of HFI data, i.e. 
bolometers are highly sensitive to any temperature signal, not only to microwaves, 
therefore any drift needs to be removed with high accuracy 

• low precision thermometers (Cernox) which serves just for housekeeping 

• Sorption cooler sensors, used to monitor the performance of the Sorption cooler 

For the rest of the chapter I will be often referring to the 4KL1 and 4KL2 sensors 
as LFI reference load sensor; it is understood that I'm referring to the 30 and 44 GHz 
reference loads, because the 70 GHz loads are located near to the top of the 4K box. 



77 




Figure 4.2.: CSL-CQM thermometers location 



Sensor 


Location 


Sampling [Hz] 


PID4N 


4K stage horns 


160 


4KL2 


Side of the 4K box 


160 


4KL1 


Side of the 4K box 


160 


Cernox_4K 


Back of the 4K box 


1/70 


tt03 


4K Cooler cold end 


160 


Di842 


LFI support 


160 


Cernox_18K 


HFI 18K plate 


1/70 


tt04 


4K Cooler pre-cooling temperature 


160 


LVHX2 


Sorption Cooler cold end on LFI 


1 


LVHX1 


Sorption Cooler cold end on HFI 


1 



4.2.4. TIME domain analysis 

The amplitude transfer function was computed in the time domain by taking the ratio 
of the peak-to-peak amplitude of the input and output fluctuations; this method has 
strong contra: 

• the transfer functions are defined on sinusoids, taking the ratio of input square 
waves amplitudes versus output triangular waves is a strong simplification 

• a slow thermal drift or high frequency noise can worsen the estimation 
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However, this method has the advantage of being very simple and fully under control, 
i.e. it is possible to check its results by eye, therefore is a good method for having a 
rough but robust estimate of the amplitude transfer functions. 



666 seconds sinusoidal input, PID Off 

This test was performed with an input square wave at LVHX1 of an amplitude of ~ 0.5K 



and a frequency of 1.5 mHz. Figures 4.3 and 4.4 show the resulting fluctuations measured 



by the high precision thermometers; figure 4.5 shows the comparison between the square 
input fluctuation at LVHX1 and LVHX2 and the outputs at the top, middle and bottom 
of the 4K box. 




0.002 



0.001 



0.000 



-0.001 



-0.002 
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Time (minutes) 
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Figure 4.3.: 4K LFI reference loads Figure 4.4.: 4K plate fluctuations (mean 
fluctuations removed) 



Table 4.1 lists the measured outputs and the computed amplitude transfer functions. 
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Figure 4.5.: 666 sec, PID OFF: comparison between LVHX, Cernox, 4K plate and 4K 
load 



80 



Sensor 


Fluctuation amplitude (K) 


Transfer function from LVHX1(K/K) 


PID4N 


0.00383 


0.00719 


PID4R 


0.00370 


0.00695 


PID4C 


0.00339 


0.00637 


4KL1 


0.00393 


0.00739 


4KL2 


0.00438 


0.00823 


LVHX2 


0.356 


0.669 


LVHX1 


0.532 


1 



Table 4.1.: Transfer function gains deduced from PACE 666sec fluctuations measurement 



100 seconds sinusoidal input, PID Off 

This test was performed with an input square wave at LVHX1 of an amplitude of ~ 0.5K 



at a frequency of 10 mHz. Figures 4.6 and 4.7 show the resulting fluctuations measured 



by the high precision thermometers; figure 4.8 shows the comparison between the square 
input fluctuation at LVHX1 and LVHX2 and the outputs at the top, middle and bottom 
of the 4K box. 
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Figure 4.6.: 4K LFI reference loads Figure 4.7.: 4K plate fluctuations (mean 
fluctuations removed) 
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Figure 4.8.: 100 sec,PID OFF: comparison between LVHX, Cernox, 4K plate and 4K 
load 
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Sensor 


Fluctuation amplitude (K) 


Transfer function from LVHX1(K/K) 


PID4N 


0.000269 


0.000384 


PID4R 


0.000294 


0.000420 


PID4C 


0.000241 


0.000345 


4KL1 


0.000336 


0.000480 


4KL2 


0.000376 


0.000537 


LVHX2 


0.690 


0.986 


LVHX1 


0.700 


1 



Table 4.2.: Transfer function gains deduced from PACE lOOsec fluctuations measurement 
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4.2.5. FREQUENCY domain analysis 



Here we discuss the method to extract the thermal transfer functions based on frequency 
domain analysis. The procedure is described step by step taking as example the 100 
seconds PID off test. 

100 sec PID off 

Selection of the stable part of the timestream The first step is to remove transients 
and drifts which could impact to the FFT. In this case the timestream is very clean and 



the drift on the thermometers is only of few mK, see figure 4.9 



In each case a trade-off was performed between removing the drift and having enough 
samples for a sufficient frequency resolution 
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Figure 4.9.: LVHX1 100 sec fluctuations 



FFT The Fast Fourier Transform module used is the FFTPACK included in Numpy 
1.0 Q the output is then divided by the sampling frequency, no other normalisation is 



1 Numpy is a multi dimensinal array and mathematical package for Python, see http: //numpy . scipy. 
org 
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Figure 4.10.: 4K LFI reference loads 
fluctuations 
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Figure 4.11.: 4K plate fluctuations 



performed because the objective is just to compute ratios. 

In the 100 seconds case the input fluctuations are sine waves while for the 666 and 
4000 case they are square waves, therefore the frequency peak is not always as clear as 



in this example, see figure 4. 12 but it is always defined. 



Computation of the transfer function The final step is to identify the peak in the 
frequency domain around the input signal frequency and take the ratio between the 
peaks. 

It was also taken into account to consider not just the peak point but also the sur- 
rounding points, by making an integral over 3 or 5 points and taking the ratio of this 
integral; however, even the nearest points are at least one order of magnitude lower, and 
already at the level of the mean of the signal, therefore they are not significant for the 
measurement; for this reason it was judged more reliable to consider just the peak point 
for taking the ratio. 



Figure 4.13 shows the peaks in a logarithmic scale. 

The results, compared to the frequency domain analyses, are displayed in table 



Sensor 


FFT 


Time domain 


Difference (%) 


PID4C 


0.000414 


0.000345 


17% 


4KL1 


0.000547 


0.000480 


12 % 



Table 4.3.: Comparison between time and frequency domain results for 100 sec PID OFF, 
the percentage difference was computed as (FFT-Time)/FFT 
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Figure 4.12.: FFT of the input LVHX1 signal and the thermometer measurement on the 
top of the 4K box, PID4C and in the middle 4KL1. The FFT is the raw 
output of the Numpy FFTPACK, no normalisation was performed 
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Figure 4.13.: FFT of the 3 signals - peaks on a logarithmic scale 
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Results 

The next table and the figures gather all the results obtained in the Fourier analysis: 



Sensor 


100s PID off 


666s PID off 


THER4KL1 


0.0005466 


0.005048 


THER4KL2 


0.0006102 


0.00562 


PID4C 


0.0004145 


0.004423 


PID4N 


0.0003987 


0.005011 


PID4R 


0.0004556 


0.004853 


THER4KH1 


0.0004585 


0.00489 


THER4KH2 


0.0004788 


0.005106 



Table 4.4.: Amplitude transfer functions, PID off 



Sensor 


100s PID on 


666s PID on 


4000 PID on 


THER4KL1 


0.0002789 


0.001344 


0.003143 


THER4KL2 


0.0003105 


0.001484 


0.003451 


PID4C 


4.898e-05 


1.512e-05 


8.703e-06 


PID4N 


5.197e-05 


1.371e-05 


5.764e-06 


PID4R 


5.499e-05 


3.204e-05 


6.475e-05 


THER4KH1 


5.3e-05 


1.225e-05 


8.208e-06 


THER4KH2 


5.609e-05 


1.843e-05 


2.583e-05 



Table 4.5.: Amplitude transfer functions, PID on 
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4.2.6. 4K timelines simulations 



The experimental results just showed were obtained with a set of PID settings still 
preliminary; the HFI team performed an optimisation of the PID settings in order to 
reduce their power consumption while still meeting the stability requirements. The 
optimisation was performed the HFI thermal software model. 

The new PID settings were then used as input to their thermal model in order to 
produce a long duration simulation (several months) of the expected fluctuations on the 
4K stage. These simulations used as input a realistic datastream of temperatures at 
LVHX1 simulated by the Sorption cooler team. 

I have used the same Fourier method for analysing this dataset; however, in this 
case I could compute the transfer function for all the range of frequencies permitted 
by the signals sampling frequencies. In fact in measurements it is necessary to give 
a strong source signal in order to disentangle effectively what component of a sensor 
measurement is due to the input source and what part instead has a different cause. 
Simulations instead are simpler and there is no doubt that all the signal we see on each 
4K sensor is actually the effect of the Sorption cooler fluctuations. 

Figure 4.14 shows the results of the comparison between the transfer functions com- 
puted from the simulations and the transfer functions from measurements. As expected 
the simulated transfer functions lie between the PID OFF and the PID ON measure- 
ments, which is the effect of the reduction of PIDs power consumption. 



4.3. Conclusion 

In this chapter I have presented the estimation of the amplitude thermal transfer func- 
tions of the HFI 4K box with respect to the Sorption cooler cold end LVHX1, with and 
without PID regulation. They have been estimated from the qualification tests held 
in CSL in 2005 using an analysis method in the Fourier domain which proved to be 
consistent with a simpler time domain method at the level of 17%. 

However, the expected behaviour in flight is different, because the setting of the PID 
thermal regulators has been optimised for reducing power consumption. The expected 
performance was simulated by the HFI team which provided a long duration simulation 
of the temperature measurements of all 4K sensors considering a realistic Sorption cooler 
input. I have then computed the expected transfer functions from the simulations: their 
comparison with the measurements is consistent; as expected, the expected transfer 
functions give a low pass filtering effect weaker than PID ON measurements but better 
then PID OFF measurements. 

The next step is the analysis of the measurements acquired during the flight model 
cryogenic testing campaign held in CSL in summer 2007. During this campaign there 
were not dedicated tests to characterise transfer functions; however, computing the am- 
plitude thermal transfer functions during a long stable acquisition test we can check 
whether the prediction of the simulations were accurate. 
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Transfer Functions comparison:HFI_95_THER_PID4N_Tl_K_C 
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Figure 4.14.: Comparison of the transfer functions computed from the simulations (black 
curve) and measurements at CQM with PID ON and PID OFF 
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4.4. Annex 
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Figure 4.16.: 666 seconds PID OFF test, complete timelines 
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5. Combined quick look analysis 



5.1. Introduction 

Quick look analysis represents the first check of instruments health and performance; 
precisely it refers to the operation of visualising the data in real-time and performing 
simple analysis tasks, like producing spectral densities, change data binning and so on. 

This analysis has been regularly performed during testing campaigns to check in real- 
time the output of the scientific channels and of housekeeping, i.e. thermometers, in- 
struments supply and control biases, power supply parameters and other. 

The HFI team used KST for this purpose, while the LFI team used Rachel just for tests 
on single RCAs and TQL (Telemetry Quick Look) for tests of the integrated instrument. 

The same tools used during the testing campaigns will be used during the commis- 
sioning phase (CPV), i.e. the first 90 days of the flight dedicated to check functionality 
and to set the instruments in the best configuration before starting the CMB survey, 
and during the rest of the mission. 

Its use during the tests, both on ground and during CPV, it fundamental in order 
to monitor in real-time what is the instrument response to the commands sent and the 
testing procedures require that the test operator checks specific scientific or housekeeping 
channels during test execution. 

During the surveys the objective of the quick look analysis is different, mainly focused 
on monitoring that no anomalies are present and that each instrument keeps its nominal 
configuration, which assures the best performance. 

Since Mission planning it was clear that also a quick look analysis involving simulta- 
neously both instruments would be very useful in order to identify possible systematic 
effects on the cross correlations; this exploits one of the big strenghts of Planck, which is 
having instruments with completely different technologies in the same focal plane; cross 
correlation can highlight issues that are not visible analysing a single instrument or can 
allow an immediate cross check on the other instrument if something unexpected is seen 
in the data stream and gives more chances of identifying its source. 

During Planck Mission, LFI and HFI scientific and housekeeping data will be ex- 
changed daily using the Exchange Format interface (EFDD), see |21j . Exchange Format 
is a specific fits file format designed to be easily transferred and managed by the other 
DPC; it is already planned that the exchange will be organised in 2 subsequent phases, 
see PQ: 

1. converted timelines will be created daily and available the day after data reception 
from the satellite, converted involves a minimal systematic effect removal and 
calibration to V for scientific channels and K for thermometers 
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2. reduced instead includes also a deeper systematic removal and a calibration of 
scientific channels to K or to W; they will be created once a week for the data of 
the previous week. 

In detail: 

• The HFI converted science timeline data is composed by the data, demodulated, 
low-pass filtered and converted to volts using the appropriate gain 

• The HFI reduced science timeline data is composed by the converted data after 
the systematic effect removal, (the headers of the files will include conversions to 
sky units for the bolometer and to Kelvin for the thermometers). 

• The LFI converted timeline is composed by differenced radiometer data after the 



best gain modulation factor (r, see section 2.3.2) application and calibrated in 
Volts. 

• The LFI reduced science timeline is composed by differenced radiometric data 
after the r factor application, after calibration to Kelvin and after the systematic 
effect removal 

Following this concept, as soon as the data of the other instrument in Exchange 
Format arrives the day after reception, an operator should manually inspect them with 
the support of a handy graphical tool. 

Moreover, this tool can be useful in later phases when a poorly understood effect for 
example on LFI data could be easily checked against HFI scientific or housekeeping data. 



5.2. KST for data crosscorrelation 
5.2.1. KST 

KST software was created by Barth Netterfield [2], an astrophysicist of the University 
of Toronto. He needed a tool to look at the data in real time during acquisition and to 
inspect quickly large datasets. 

The software is coded in C++ and based on the QT graphical libraries and the KDE 
libraries, it was released as Free Software under GPL and is now released together with 
the K Desktop Environment. 

It is a GUI oriented plotting and analysis tool focused on working with large datasets or 
directly interfaced with streaming acquisition electronics; the implementation is focused 
on the GUI being quick and easy to use. It also supports plugins in order to extend its 
features. 

It was since then used by the Cosmic Microwave detection balloon experiment Boomerang 
([H])p4j) as a tool for graphical plotting and navigating the data, see an example in 
figure |5.1| 

The Canadian Space Agency funded the project in order to support the ESA Planck 
Mission with 2 main objectives: 
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Figure 5.1.: Screenshot of KST plotting BOOMERANG data 
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1. Being the HFI quick look analysis software used for plotting and inspecting live the 
data coming from the satellite and the tool for quickly inspecting data afterwards. 

2. Being the designated tool for data cross correlation between HFI and LFI to be 
used to manage data in Exchange Format 

5.2.2. Concept of quick look analysis sessions 

KST is very efficient for loading data and producing plots, however manually loading 
data, adjusting legends and labels is quite tedious to be performed routinely. 

Therefore the need of automating this tasks has been fulfilled by preparing a config- 
urable set of scripts which will manage the repetitive part of data visualisation. 

Beside this, having predefined analysis sessions helps standardising the way different 
operators look at the data, improving robustness and reliability of data analysis. 

Another interesting KST feature is that once a part of a dataset, for example two 
hours of data, including amplitude spectral densities or histograms, it is visualised, it 
is easy to switch to the next hours back and forth just by pressing buttons; this allows 
to set up even a complex analysis environment and then apply it sequentially to the 
complete duration of the dataset. 

5.2.3. Implementation 

KST supports a scripting language and an interactive console which is available at the 
bottom of KST GUI; the implementation is based on javascript and binds most of the 
C++ objects which form the KST window. 

It is possible to run command line functions interactively on KSTscript console, and 
prepare scripts offline that can be run and take arguments from this console. 

KSTscript allows to load data, create windows, make plots, modify labels, create 
power spectra and so on. 

The first objective was to create an abstraction layer for accessing the data, in order 
to be able to build analysis sessions without worrying about how the data are stored. 

This abstraction layer is a javascript library named PlanckData I wrote using an object 
oriented structure, it allows the access to different data format with the same syntax. 

For now the available dataset interfaces are: 

• Level S 24 hour simulation: HFI and LFI simulated timestreams produced by the 
Planck simulator, called Level S 

• Exchange format: ability to read fits files in Exchange format structured in prede- 
fined folders 

It is important to note that data access itself is already built in into KST, what this 
library allows is to call the data from javascript with a very simple syntax, typically, 
once defined what data source we are using, the syntax is: 
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source = new PlanckData (" 24 hLevelS " ) ; 

time = source . ch [" LFI23a "]. getTime ( start , stop , binning ) ; 
signal = source . ch [" LFI23a "]. getSig ( start , stop , binning ) ; 

This example is based on the 24hlevelS dataset which makes use of the old LFI naming 
conventions; this dataset contains data already differenced, therefore each RCA has an 
a and a b streams. 

The library also provides consistent labelling to all sessions, sampling frequencies 
computations and data binning. 

Scripts do not directly manage the data, they just control the behaviour of the C++ 
objects which do the hard work. 

5.2.4. Base sessions 

Base sessions perform atomic actions on KST, and they are the basis for building ad- 
vanced sessions. 

Each base session has an inline help and a test function in order to improve user 
friendliness. 

Sessions already implemented includes data loading and plotting functions for: 

• data streams 

• amplitude spectral densities 

• histograms 

• simple correlation 

• correlation function 

• autocorrelation 

• cross correlation 

• cross power spectrum 

User interface is managed in 3 different ways: 

• command line interface, i.e. call a function with arguments: 

tod({ch:[ '2700 ', '2410 ', '1911 ', '2201 '] , t0:0, dt:2, bin:0, over:0}); 

where ch is channel list, tO begin time, dt time interval, bin binning and over 
defines whether each curve should be on its own plot or overplotted. 

• a pop up dialogue interface, which asks the user the needed arguments with hints 

• a complete GUI, with a table of all available channels to be selected by clicking on 
tickboxes 
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tad({ch : ri^-^bv 1^1-273 v jag- 1*' ,'317-3'] 
kst* 



to : Q, d1 : 1. bin:[J, a*er:L>); 



Figure 5.2.: Overplotting of HFI and LFI simulated data streams produced with tod 
session, in the console there is the function call 



Basic sessions have been already developed and tested on Level S test data and on an 
Exchange format dataset and are fully functional; they are suitable for supporting the 
activity of an operator looking at the data, but they are not proper analysis sessions by 
themselves. 



5.2.5. Advanced sessions 

The role of advanced sessions is different, their aim is to be a good, even if partial, view 
of a possible specific issue; therefore they should be self standing and should contain all 
the information needed for understanding if an issue is present or not. 

This type of session could also be used as a starting point for reporting, i.e. the 
operator, after having looked through the plots, could log that there is no evidence of 
unexpected signatures and mark that the data range is clean. 

In this context I'm evaluating the possibility of exploiting the available interface with 
the ELOG software, which is an electronic logbook which was used by HFI during test 
campaigns to log all the operations performed on the instrument directly from within 
KST, being also able to save a .kst file or a KST screenshot. 

This could be handy to keep track of all the analysis performed on data cross corre- 
lation and having the possibility to easily search through the electronic logbook to go 
back to a previously identified issue. 
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Figure 5.3.: All possible crosscorrelation function plots of the previously selected chan- 
nels produced by xcorr 
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Figure 5.4.: Overplotted amplitude spectral densities produced by asd 
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5.2.6. List of proposed advanced analysis sessions 

Before discussing HFI and LFI data cross correlation, we first discuss simpler sessions, 
based on data from a single instrument. This will help understanding the sessions 
structure and the main feature present in the data. 

LFI-only sessions 

• Power group based: 

LFI is biased thanks to four different powergroups; one of the first analysis tasks 
is to check if there are common effects on the same powergroup. Here follows a 
list showing the correspondence between the RCAs and the powergroups: 

- powergroup 1: 18,26 

- powergroup 2: 19,20,28 

- powergroup 3: 21,22,24,27 

- powergroup 4: 23,25 

In order to simplify the analysis it would be better to choose the same frequency, 
therefore 70 GHz, whose RCA are spread on all powergroups. 

Proposed sessions: 

1. 18,19,21,23: crosscorrelation between different powergroups 

2. 19,20: crosscorrelation on the same powergroup 

3. 21,22: crosscorrelation on the same powergroup 

Total of 8 plots, which could be accommodated on a single KST window, we can 
choose to plot just one channel for each RCA listed using differenced data 

• Frequency based: 

26,23,28 different frequencies and different powergroups in order not to see effects 
in the same powergroup 

3 plots 

• RCA internal correlation 

The correlation between the channels of the same RCA should show the strong 
correlation of the channels of the same radiometer due to the 1// amplifiers gain 
fluctuations; while the correlation with the channels of the other radiometer is just 
due to the sky signal. The number of plots for checking all the correlations for 
each RCA are: 

- 6 single channel; 

- 1 combined diodes. 
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HFI-only sessions 

Considering a channel for frequency and type (i.e. spider web or polarisation sensitive 
bolometers), we have 9 channels: 



1. 


100- 


lb 


2. 


143- 


la 


3. 


143- 


5 


4. 


217- 


5a 
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217- 


2 


6. 


353- 


3a 


7. 


353- 


7 


8. 


545- 


3 


9. 


857- 


3 



HFI channels readout electronics is organised in groups of six channels; the channels 
on the same group share the same FPGA (Field-programmable gate array) controls the 
channel settings and performs the AC demodulation. 

A possible consistent choice of sessions, which was discussed with HFI experts, can 
be: 

• 2 channels at different frequency on the same belt, in order to identify possible 
electrical effect on the belts 

• 2 nearby channels in order to check for possible correlations due to thermal coupling 

• both bolometers of a Polarisation Sensitive Bolometer 

• a channel for each frequency 

The objective of the KST sessions is to be quickly analysed by a human operator, 
therefore it is necessary to choose among these possible cross correlations a consistent 
subset with no more than two windows and less than nine plots each. 



HFI-LFI cross correlation sessions 



Satellite power supply monitoring 

A very simple correlation analysis involves monitoring the electrical power supply 
line connecting the scientific instruments to the satellite. 

In case of anomalies in the instrument outputs, it could be useful to have a session 
ready to cross correlate voltage and current measurements of the supply line. 



Tables 5.1 and 5.2 summarise the HFI and LFI housekeeping parameters that can 



be used in this correlation analysis. 
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Table 5.1.: HFI supply line parameters 

Parameter Code Description 

HM013250 DPU + 5V (V) 

HM014250 DPU - 15V (V) 

HM015250 DPU +15V (V) 

HM016250 DPU +2.5V (V) 

HM006250 DPU Intensity of +5V (mA) 

HM007250 DPU Intensity of -15V (mA) 

HM008250 DPU Intensity of +15V (mA) 

HM009250 DPU Intensity of +25V (mA) 



Parameter Code 



Table 5.2.: LFI supply line parameters 
Description 



LM214342 
LM215342 
LM216342 
LM218342 
LM219342 
LM220342 



VCC Current (A) 
+15V Current (mA) 
-15V Current (mA) 
VCC Voltage (V) 
+15V Voltage (V) 
-15V Voltage (V) 



• Crosscorrelation between 70 & 100 GHz data 

This analysis can be performed by correlating data from two LFI 70 GHz and two 
HFI 100 GHz horns, on symmetric positions on the focal plane; cross correlating 
for example a 70 GHz channel with a 100 GHz very near and another on the 
opposite side of the focal plane can highlight correlations due to their proximity, 
e.g. something related to cross talk or to non idealities in a region of the telescope 
reflectors. 

All the possible crosscorrelations are just six and can fit in a single KST window 
or being overplotted. 

• Large frequency coverage cross correlation 

A cross correlation involving a wider set of Planck channels could check for signa- 
tures of unexpected correlated effects coming, for example, from the telescope. 

In order not to increase too much the time the operator needs to go through this 
session, a subset of channels per frequency can be: 

— a LFI 44 GHz channel differenced 
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- a LFI 70 GHz channel differenced 

- a HFI 100 GHz channel 

- a HFI 143 GHz channel 

- a HFI 217 GHz channel 

All the possible cross correlations, which could be either performed using the cor- 
relation function or using cross spectra, use 10 plots. 

• Polarisation dedicated sessions 

Possible polarisation dependent effects could be investigated by correlating the 
polarisation measurement from LFI 70 GHz and HFI 100 GHz. 

For LFI it is necessary to consider two radiometers connected to the orthogonal 
polarisations of the same horn, for each of them taking the difference of sky versus 
load and calibrate to sky units; at that point polarisation is measured by taking 
the difference of the two streams. This procedure gives a first estimation of the 
polarised signal, for a deeper analysis it is necessary to combine the data of the 
two horns symmetric in the focal plane. 

The corresponding procedure should be performed on HFI 100GHz, exploiting the 
data of the orthogonal Polarisation Sensitive Bolometers. 

The reduced exchange data are calibrated to Kelvin units, LFI data are already 
differenced, therefore this analysis can give a rough estimation of the polarised 
signal as a function of time. 

The previous sessions mainly examined issues that shouldn't be present in nominal 
conditions; the only correlations should come from the sky signal, which will be removed 
by a dedicated algorithm already implemented into the data analysis pipeline. 

The next session instead is fully dedicated to thermal issues that are already expected 
to be present; therefore the aim is to check that they comply with requirements and 
identify any anomaly. 

• sessions dedicated to the 4K stage 

LFI reference loads are attached to the HFI 4K shield; correlations between the LFI 
scientific signal and 4K loads fluctuations are already top priority for correcting 
the observation results: it is already planned that the data analysis pipeline will 
ingest the 4K thermometers data in order to clean the LFI scientific signal from 
these fluctuations. 

The objective of an analysis with KST instead is focused on involving also HFI 
scientific channels, because the focal plane, where feed horns are located, is part 
of the 4K stage, and its fluctuations could impact the bolometric output. 

Thermal fluctuations are very different depending on the position on the 4K shell; 
the interface to the 4K cooler is at the bottom of the structure, where the fluctu- 
ations are stronger; on the top plate, which is the HFI focal plane, instead, there 
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is a strict thermal regulation performed by PID controlled heaters, therefore the 
temperature is very stable; in the middle, where LFI 30 and 44 GHz reference 
loads are located, there is an intermediate situation which is difficult to model due 
to the complex design. 

Therefore there are two possible scenarios: 

— top thermal plate fluctuations: in this case the it would be useful to correlate 
a 70 GHz channel differenced data stream, an HFI thermal sensor located on 
the top plate and a 100 GHz channel. 

— LFI 30/44 GHz load fluctuations: in this case a 30 GHz channel, a 44 GHz 
one, a 100 GHz and one of the HFI 4K belt thermometers (4KL1/2) should 
be used. 



5.3. Preliminary tests on real data 

Planck satellite was tested in Liege Space Centre (CSL) during Summer 2008; this was 
the first time HFI and LFI worked together integrated into the satellite and in flight-like 
cryogenic conditions. 

It was not in the test plans to produce exchange format and perform cross correlations 
between the two instruments data, 

but the long duration test, where HFI and LFI acquired data in stable conditions 
for 24h was a really good opportunity for testing both exchange format production and 
cross correlation analysis. 

Unfortunately, due to bugs in the exchange format conversion software, exchange 
format data were made available after the conclusion of CSL campaign, after much data 
were already been analysed with other tools. 

However it is very interesting to use this dataset for testing the concept of sessions 
and also further analyse this dataset which still shows few unclear effects. 

The main problem we are facing is that the exchange format data available now are 
just the converted streams, which means they are calibrated to volts, while most of the 
analysis sessions described in this chapter relies on scientific data calibrated in Kelvin 
units. This is the only way we can really compare data from different channels. 

In the next future probably CSL data will be made available also as reduced timestreams, 
until then we can start testing the environment and few sessions. 



In figure [53] HFI 100 GHz lb is overplotted to LFI 19 70 GHz S-ll, HFI bolometers 
show a strong drift, probably due to a thermal drift in the cold HFI stages. Therefore 
I have subtracted this drift by differencing the HFI 100 GHz lb channel with one of 
the dark bolometers, which are bolometers equivalent to the others but that are not 
connected to any horn and are therefore very useful to distinguish between thermal 
effects and sky signal. 

It shows a few hours period when the thermal shield was refilled with helium, which 
caused is a distinct temperature change recorded by all HFI and LFI channels. Then 
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Figure 5.5.: LFI19S-11 and HFI 100-1B channels overplotted during CSL long duration 
test 

I have removed the mean and applied an arbitrary scale factor in order to mimic the 
effect of a calibration to Kelvin and obtain a similar step height. 

Figure [576] shows the cross correlation function between the data from these channels; 
as expected there is a very strong correlation around offset, which means that the 
signal is very synchronous on both channels. 

5.4. Conclusion 

This project kst4planck, developed by Rodrigo Leonardi (UCSB), Luis Mendes (ESA), 
Andrew Walker (UBC) and myself, is focused on delivering a set of predefined quick 
analysis sessions in order to simplify the process of cross correlating HFI and LFI data; 
the objective is to support the removal of systematics and to look for unexpected issues 
during the survey. 

At this point there is already a set of working base sessions, its development showed 
that the session concept is really applicable to KST and technically successful. 

The next important step is to demonstrate its validity also on the field, by concentrat- 
ing on the CSL long duration test and improve the knowledge about the environmental 
conditions by exploiting the crosscorrelation of scientific data and housekeeping of both 
instruments. 

This is a fundamental step on the path for delivering a flight-ready tool both to the 
LFI and the HFI DPCs; kst4planck will be available among other tools and will be 
specifically dedicated (but not only) to exchange format data crosscorrelation. 
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Figure 5.6.: Crosscorrelation of LFI19S-11 and HFI 100-1B channels 
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6. Conclusion 



The work I have discussed in this PhD thesis is focused on three main topics: LFI 
radiometers software modelling, 4K stage thermal transfer functions and HFI-LFI cross- 
correlation. 

This work has been performed in the context of the ESA Planck Mission, which has the 
objective of providing the definitive measurement of the Cosmic Microwave anisotropies 
temperature down to the smallest scales, a good measurement of anisotropies E-mode 
polarisation and possibly the first detection of B-mode polarisation. 

QIMP model and radiometers bandpasses During my PhD I have developed QUCS 
Integrator of Measured Performances (QIMP), a software model which relies on measured 
frequency response of the LFI radiometers component and on an analytical waveguide 
simulator in order to simulate the LFI bandpasses channel by channel. 

Precise knowledge of radiometer bandpasses is fundamental for effective foregrounds 
removal. Foregrounds removal algorithms work on maps at different frequencies by 
exploiting the different spectral indexes of the foreground components. However, each 
HFI and LFI channel has a broadband response, and each map is the convolution of 
the spectra of foregrounds and CMB with the instruments frequency response. Thus, 
an iterative process involving component separation and maps rescanning with realistic 
bandpasses is needed to take this effect into account. Moreover, detailed bandpasses 
mismatch characterisation of the two radiometers of the same horn allows the study of 
the leakage of foregrounds temperature signal into polarisation signal, due to the fact 
that the calibration is tailored to the spectrum of the CMB, while foregrounds have 
different spectra. 

LFI bandpasses were measured on dedicated tests during the integrated RCAs cam- 
paigns; however a part of these measurements showed strong systematic effects due 
to the test setup (30 GHz and 70 GHz) or had an insufficient frequency coverage (30 
GHz). The QIMP model can compensate for this issue by providing bandpass estima- 
tions based on frequency response tests on the single RCA components and an analytic 
waveguide simulator. QIMP results were validated against measurements and showed a 
good agreement, with the exception of a discrepancy on 44 GHz channels which is still 
under investigation. 

Single components performance tests (which are less prone to systematic effects be- 
cause of simpler setups) and QIMP simulations provide the best knowledge of LFI ra- 
diometers frequency response. 

4K stage thermal transfer functions I have analysed the thermal transfer functions 
between the Sorption cooler cold end LVHX1 and all thermometers on the 4K stage, see 
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chapter |4j 

I had the opportunity of computing the transfer functions on data acquired during 
the HFI CQM cryogenic tests in Liege Space Centre with and without PID thermal 
regulation and compare these results with the thermal functions computed from the 
simulations performed with the HFI thermal model after PID settings optimisation; the 
results are consistent considering that the optimisation aimed at reducing PIDs power 
consumption. 

4K stage thermal transfer functions are fundamental for LFI because they allows to 
compute what is the effect of Sorption cooler and 4K cooler fluctuations on the 4K 
reference loads. The requirement on the spectrum of the 4K load thermometers is 
10]xK/ y^[Hz) for frequencies over WmHz and W0[iK/y/(Hz) on lower frequencies; the 
last estimation of this effect, dated June 2006, showed that for 30 and 44 GHz RCAs the 
requisite was not met, while for 70 GHz it was met. In the next future the fluctuations 
will be analysed both by using the latest transfer functions and by directly computing 
the spectrum of the thermometers, in order to check what is the impact on maps of these 
fluctuations. 

HFI- LFI cross-correlation Starting from September 2007 I have led the development 
of the kst4planck project, with the objective of designing and implementing quick look 
analysis sessions for the cross-correlation of HFI and LFI scientific and housekeeping 
data in Exchange format, see chapter [5] 

Analysis sessions are implemented in KST, a data visualisation software; we have 
developed and tested the basic sessions for loading, displaying and performing simple 
analysis, like spectra, cross-correlation, cross power spectra. 

Basic sessions are the elements for building high level sessions, i.e. a integrated set 
of windows and plots focusing on a specific issue, that will allow an operator to quickly 
navigate the data for identifying unexpected systematic effects in the cross-correlations. 
Among the huge amount of data and possible cross-correlations we have identified a set 
of useful cross-correlation sessions that we are now ready to implement. 

HFI-LFI cross-correlation is a topic of growing interest in the Planck community; 
Planck peculiarity of having two very different detection technologies, and therefore 
different systematic effects and different signatures in the data of external systematics, 
is a strong opportunity for identifying systematic effects and remove their effect from 
the scientific data. 

6.1. Future work 

QIMP QIMP still has strong margins of improvement: first of all thanks to the dedi- 
cated large band and section by section bandpass tests which were performed by Stefania 
De Nardo in her master thesis (|12j) on a 30 GHz spare RCA. A detailed comparison 
between the test results and the model simulations will improve the model bandpass 
characterisation; the improved model will then be used to run again simulations on the 
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flight bandpasses. Moreover, 44 GHz test data will be analysed again in order to figure 
out if this could have caused the bandpass discrepancy. 

QIMP can also be used to model other instrument performances, its use is very inter- 
esting because we can have a deeper insight on what component inside the RCA is the 
reason of a specific behaviour. Follows a list of interesting analysis topics that can be 
covered by QIMP: 

• implement the non linear model for 30 and 44 GHz RCA by including into the 
BEM model the compression factor measured during the test campaigns 

• use the results of the measurements of the radiometer output as a function of input 
temperature (LIS) in order to validate the absolute output of the 30 and 44 GHz 
radiometers 

• estimate LFI radiative emission: an analysis on LFI radiative emission was per- 
formed by Cristian Franceschet in his master thesis [TE], he used a radiometer 
model based on LFI nominal performance and showed that LFI radiative emission 
towards HFI is about two orders of magnitude below requirements; it would be 
interesting to analyse now LFI emission with the real radiometers performance 

• improve the modelling of the bandpass response for the reference load including 
second order effects like spillover, return loss and eventually straylight (in collab- 
oration with Francesco Cuttaia) 

Thermal transfer functions The next important step for improving the knowledge of 
the thermal behaviour of the 4K stage is the estimation of the thermal transfer functions 
in the integrated flight model satellite cryogenic tests performed at the Liege Space 
Centre in 2008. Thanks to this we can check the consistency with the transfer functions 
predicted by the simulations. The most suitable test for this analysis is the long duration 
test, where the instruments were left in stable nominal acquisition mode for more than 
24 hours. 

For concluding this study it would be interesting to fit the transfer function with a 
simple low pass filter in order to have a rough estimation of the thermal time constants 
of the system and to have a simple representation of the thermal behaviour of the 4K 
shell, that could be used in place of the full response when a quick estimation is needed. 

HFI-LFI cross-correlation kst4planck is still in a early stage of development, the first 
task is the implementation of the high level sessions, but this is just the beginning of a 
deep test work based on the CSL long duration test data aimed at refining the design of 
the high level sessions in order to be ready for the operations. 

The main focus, however, is not on the tool development, but on training people on 
understanding what we can learn by cross-correlation, where to search for systematics 
and how to better exploit this combined analysis. In fact, KST is just one of the tools for 
engaging cross-correlation, in the near future it will be important to widen the point of 
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view and understand how LIFE (i.e. the LFI radiometers data analysis software), KST 
and the DPC pipeline can be used to take the most out of HFI-LFI cross-correlation. 
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A. Scattering parameters 



A very useful technique for modelling radiometric components is through an S-parameters 
2-ports component. 

S-parameters, or Scattering parameters relate the amplitude and phase of reflected 
and output signal to the input signal as a function of frequency. In order to characterise 
the frequency response of a 2-ports component, four S-parameters need to be measured 
all over the bandwidth: 

511 reflected wave amplitude(port 1) /input wave amplitude(port 1) 

521 output wave(port 2) /input wave(into port 1) 

512 output wave(port 1) /input wave(into port 2) 

522 reflected wave(port 2) /input wave(port 2) 

Therefore the relation between input and output voltages can be expressed in matrix 
notation, see 1 30 1 : 



yout 
yout 



5*21 



£12 
S22 



yin 



(A.l) 



Return and insertion losses 

relation is: 



Return loss (RL) is sometimes used instead of Su, their 



RL = 20log w (S u )dB 



(A.2) 



Insertion loss is instead very different from S21, Insertion loss refers just to the losses 
due to internal dissipation, while S21 is the ratio between output and input amplitude, 
therefore it includes both the effect of reflection and the effect of internal dissipation. 
Their relation is: 



IL 



-10 log 



IS' 



10 



1 



21 1 
Pi] 



(A.3) 
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B. Microwave Waveguides analytical 
treatment 



Electromagnetic wave guides are hollow cavity enclosed by metallic walls. The propa- 
gation of electromagnetic waves can be calculated by solving the Maxwell equations in 
the cavity: 

V x ^ = -jufjt (B.l) 

V x ^ = jujet (B.2) 

where E is the electric field, the magnetic field, uj the angular frequency, /i is the 
magnetic permeability, e is electrical permittivity. These equations have multiple solu- 
tions, or modes, which are eigenfunctions of the equation system. Each mode is therefore 
characterised by an eigenvalue, which corresponds to the axial propagation velocity of 
the wave in the guide. 



Let us consider the special case of rectangular waveguides (see figure B.l). Because 
this waveguide is invariant under translation along z axis, the variation of the electric 
field can be written as E(x, y, z) = E(x, y) [Ae' yz + Be* 2 ] , where 7 is the propagation 
constant, and A, B two arbitrary complex constants. This waveguide structure allows 
propagation of two kinds of modes [30] : 

• TE m „ mode (Transverse Electric) which have no electric field in the direction of 
propagation. 

• TM mn modes (Transverse Magnetic) which have no magnetic field in the direction 
of propagation. 

Moreover it can be shown that modes exist only beyond a frequency called cut-off fre- 
quency. Below this frequency propagation of this mode cannot occur apart from evanes- 
cent modes that decay exponentially. For a rectangular waveguide this frequency is 
defined by: 

c fmir\ 2 /n7r\ 2 . . 



where c is the speed of light in vacuum. However, we will limit to the computation 
of the fundamental mode propagation, which is the only one necessary to characterise 
Planck LFI waveguides. Moreover, higher modes need the knowledge of excitation and 
are usually simulated using a full wave simulation (like FEM simulation or method of 
moments) . 
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Figure B.I.: Rectangular wave guide schematic 



The propagation constant 7 n ,m can be decomposed in two terms one propagation term 
called (3 and one loss term called a: 

7™ = j 'mn OLmn (-^-4) 

Using the free space wavenumber k = uj^/JIe and k Cmn = ^r^/£pf Cn m-> ^ could be show 
that: 

Pnm = yjtf - k 2 Cnm (B.5) 

Loss term can be decomposed into a dielectric component and a resistive component 
(a = a c + ad)- The dielectric loss can be computed for all modes and has the following 
form: 

k 2 tan 5 , 
a d = (B.6) 

where 5 is the material loss factor. 

The resistive loss, instead, is mode-dependent. For TEio mode we have: 

R s is the metal surface resistance which is linked to the metal conductivity, p, by: 

R s = y^JJTp (B.8) 
Another important parameter to compute is the waveguide impedance: 

Zte = j (B.9) 
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Where 77 is the vacuum impedance (rj ~ 1207T ~ 3770) 

Using this model a transmission line model (or S-parameters model) can be delined, 
where the S matrix is give by: 

ofS= [ ; o»-4>-**> (B - io) 

» n n ' 

where I is waveguide length and n is given by: 

n = 2cosh(7Z) + ( z + 1 J sinh(7/) (B.ll) 



Z c is an arbitrary chosen global circuit parameter called normalisation impedance. For 
ease of computation, it is usually taken for simulation purpose as Zq = I fl. 
Noise is computed easily from S-parameters using Bosma's theorem |39|. 
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C. Plots of all the comparisons between 
measured and QIMP bandpasses 



The following plots show the comparison between the bandpasses directly measured with 
a Swept source test and the simulated bandpasses obtained by QIMP integrating the 
frequency response measurements performed on the single components. 

On abscissa the frequency is expressed in GHz units, while in ordinate the bandapass 
consists of the ratio between output voltage in Volts and input power in Watts in the 
logarithmic units: 10 log 10([V]/[W]). 

For 30 and 44 GHz channels, the bandpass measurements are just relative to the 
central frequency, therefore their integral was normalized to the integral of the model. 



C.l. 70 GHz 

C.l.l. RCA 18 
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C.1.2. RCA 19 
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c. 



1.3. RCA 20 
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C.1.4. RCA 21 
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C.1.5. RCA 22 
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C.1.6. RCA 23 
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C.2. 44 GHz 

C.2.1. RCA 24 
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Figure C.7.: RCA 24 comparison between measurements (black) and model(blue) results 
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C.2.2. RCA 25 
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C.2.3. RCA 26 
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C.3. 30 GHz 

C.3.1. RCA 27 
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Figure C.10.: RCA 27 comparison between measurements (black) and model(blue) 
results 
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C.3.2. RCA 28 
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Figure C.H.: RCA 28 comparison between measurements (black) and model (blue) 
results 
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